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is  Phase  I  of  a  three  phase  progran. 

This  report  (Volume  IV)  presents  the  accomplishments  of  the  task  en¬ 
titled  "Initial  Quality  Representation"  of  this  program.  Four  other  vol¬ 
umes  are  written  to  describe  the  summary  of  and  the  progress  made  in 
Phase  I.  They  are: 

Volume  I  -  Phase  I  Summary 

Volume  II  -  Durability  Analysis:  State-of-the-Art  Assessment 

.Volume  III  -  Structural  Durability  Survey:  State-of-the-Art 
Assessment 

Volume  V  -  Durability  Analysis  Methodology  Development 

This  report  is  published  only  for  the  exchange  and  stimulation  of 
ideas.  As  such,  the  views  expressed  herein  are  not  necessarily  those  of 
the  United  States  Air  Force  or  Air  Force  Flight  Dynamics  Laboratory. 
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SECTION  I 


INTRODUCTION 

The  initial  fatigue  quality  of  the  durability  critical 
parts  of  an  aircraft  structure  is  one  of  the  key  input 
parameters  for  the  durability  analysis  to  be  performed  in 
the  present  investigation.  Indeed,  in  the  present  investi¬ 
gation,  we  consider  that  either  the  EIFS  (Equivalent  Initial 
Flaw  Size)  or  TTCI  (Time  To  Crack  Initiaion)  represents 
such  an  initial  quality  that  controls  the  durability  of  air¬ 
craft.  A  statistical  characterization  of  the  EIFS  data  is 
attempted  by  fitting  various  probability  distribution  func¬ 
tions  thereto.  The  EIFS  data  considered  here  are  from  the 
"Fastener  Hole  Quality"  program  [1]  and  are  listed  in 
Table  1.  These  data  consist  of  data  sets  XQPF,  XWPF  and 
WPF  of  size  37,  37  and  38,  respectively. 

To  be  more  specific,  these  data  sets  are  examined  to 
see  if  they  fit  (a)  the  distributions  of  the  Johnson  family 
including  the  long-normal  distribution,  (b)  the  Weibull  dis¬ 
tribution,  (c)  the  distributions  of  the  Pearson  family, 

(d)  the  asymptotic  distributions  of  largest  values  and 

(e)  the  TTCI  compatible  distributions.  It  is  found  that, 
with  respect  to  the  data  examined,  the  Weibull  compatible, 
log-normal  compatible,  and  second  asymptotic  distributions 
provide  the  best  fit. 


1 


Table  1  Ordered  Observations,  Their  Plotting  Positions, 
and  Some  Sample  Statistics  (EIFS  In  mils) 


TEST  SERIES  *  XQPF  TEST  SERIES  =  XWPF  .  TEST  SERIES  =  WPF 


NO, 

E IE  S  CUM 

.PROB. 

NO. 

EIFS  CUM 

.PROB. 

NO. 

EIFS  CUM 

.PROB. 

l 

.026 

.0263 

1 

.093 

.0263 

1 

.190 

.0256 

2 

.031 

.0526 

2 

.195 

.0526 

2 

.236 

.0513 

3 

.097 

.0789 

3 

.195 

.0789 

3 

.280 

.0769 

4 

.050 

.1053 

9 

.160 

.1053 

9 

.280 

.1026 

5 

.058 

.1316 

5 

.175 

.1316 

5 

.320 

.1262 

6 

.059 

.1579 

6 

.180 

.1579 

6 

.367 

.1538 

7 

.060 

.1892 

7 

.160 

.1892 

7 

.367 

.1795 

B 

.063 

.2105 

8 

.190 

.2105 

8 

.920 

.2051 

9 

.063 

.2368 

9 

.190 

.2368 

9 

.950 

.2308 

10 

.  068 

.2632 

10 

.210 

.2632 

10 

.950 

.2569 

11 

.066 

.2895 

11 

.210 

.2695 

11 

.950 

.2821 

12 

.090 

.3158 

12 

.290 

.3158 

12 

.982 

.3077 

13 

.096 

.3921 

13 

.290 

.3921 

13 

.518 

.3333 

19 

.105 

.3689 

19 

.290 

.3689 

19 

.520 

.3590 

lb 

.110 

.3997 

15 

.290 

.3997 

15 

.520 

.3896 

16 

.113 

.9211 

16 

.250 

.9211 

16 

.560 

.9103 

17 

.120 

.99  79 

17 

.290 

.9979 

17 

.560 

.9359 

IB 

.130 

.9737 

18 

.295 

.9737 

18 

.600 

.9615 

19 

.  150 

.5000 

19 

.295 

.5000 

19 

.697 

.9872 

20 

.170 

.5263 

20 

.295 

.5263 

20 

.697 

.5128 

21 

.200 

.5526 

21 

.295 

.5526 

21 

.650 

.5385 

22 

.220 

.5769 

22 

.330 

.5789 

22 

.698 

.5691 

23 

.270 

.6053 

23 

.330 

.6053 

23 

.698 

.5897 

29 

.300 

.6316 

29 

.330 

.6316 

29 

.698 

.6159 

25 

.536 

.6579 

25 

.920 

.6579 

25 

.698 

.6910 

26 

.536 

.6892 

26 

.920 

.6892 

26 

.698 

.6667 

27 

.611 

.7105 

27 

.920 

.7105 

27 

.698 

.6923 

28 

.612 

.7368 

26 

.970 

.7368 

28 

.759 

.7179 

29 

.650 

.7632 

29 

.9/0 

.7632 

29 

.759 

.7936 

30 

1.090 

.7895 

30 

.590 

.7895 

30 

.817 

.7692 

31 

1.090 

.8158 

31 

.612 

.8158 

31 

1.090 

•  7949 

32 

1.090 

.8921 

32 

.700 

.8921 

32 

1.190 

.8205 

33 

1.100 

.6689 

33 

.810 

.8689 

33 

1.250 

.8962 

39 

1.190 

.8997 

39 

.610 

.8997 

39 

1.250 

.8718 

35 

1 .290 

.9211 

35 

.870 

.9211 

35 

1.990 

.89/4 

36 

3. '0  00 

.9979 

36 

.990 

.9979 

36 

1  .690 

.9231 

37 

7.700 

.9737 

37 

1.280 

.9737 

37 

2.730 

.9487 

38 

3.830 

.9744 

MEANS  .6233 

MEANS 

.3868 

MEANS 

.7986 

AM2s  1.7214 

AM2  = 

.0694 

AM2= 

.4660 

AM3=  9.8704 

AM3= 

.0285 

AM3= 

.9179 

AM4S68. 7061 

AM4  = 

.0241 

AM4  = 

2.6276 

STDVs  1.3120 

ST0V  = 

.2634 

STDV  = 

.6826 

SUR1 (81 )s  4.3701 

SORT (81)= 

1 

.5598 

SORT (B  1 )  = 

2.8860 

62=23.1853 

62  = 

5 

.0060 

82  = 

12.1022 

61=19.0982 

61  = 

2 

.4331 

61  = 

8.3292 
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S  E"C  T  I  0  N  II 


STATISTICAL  ANALYSIS  PROCEDURES 

2.1  Standard  Measure  of  Skewness  ^B^  and  Standardized 
Measure  of  Peakedness  (Kurtosis)  B2 

With  respect  to  a  random  variable  X,  write  y^  for  the  i-th  moment  about 
the  origin  and  y.  for  the  i-th  moment  about  the  mean; 

yj  =  E{X7},  y.  =  E{(X  -  E[X])7}  (1) 

The  standardized  measures  of  skewness  /B^  and  of  peakedness  B2  are  then  91*" 
ven  by 

^  =  y3/(4Ip3  B2  ■  y4/(y2)2  (2) 

The  estimations  of  the  moments  can  be  made  on  the  basis  of  the  observed  da¬ 
ta  consisting  of  a  sample  size  n;  x^,  x2 . xn.  In  fact,  introducing 

the  sample  moments  ml  and  m..  as 

ml  =  (  l  x7)/n  m.  =  I  (x.  -  mi)7/n  (3) 

1  k=l  *  1  k=l  K  7 

the  population  moments  yl  and  y..  may  be  estimated  by  their  corresponding 
sample  moments  mj  and  m^  and  also  skewness  /B^  and  peakedness  B2  respec¬ 
tively  by  and  b2: 

7b^  =  m3/(v€^)3  b2  =  m4/(m2)2  (4) 

2.2  Use  of  the  B^-B2  Plane  for  the 
Selection  of  Distribution  Functions 

A  large  number  of  distribution  functions  may  be  considered  for  the  cur¬ 
rent  investigation  to  characterize  the  EIFS  data.  An  engineering  approach 
has  been  used  here  to  select  from  these  distribution  functions  only  those 
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distributions  that  have  the  prospect  of  passing  further  tests  of  goodness- 
of-fit,  while  eliminating  those  that  are  obviously  incompatible  with  the 
given  EIFS  data.  The  approach  is  to  use  the  plane  as  shown  in  Fig. 

1  [2]  where  each  of  these  distributions  may  be  identified  either  as  a 
point,  curve,  or  region.  Although  Fig.  1  indicates  the  rela¬ 
tionships  between  and  some  the  better  known  distribution 

functions,  they  have  been  established  analytically  and  do  not  necessarily 
Indicate  the  relationships  between  the  estimates  b^  and  b^  of  and 
Nevertheless,  Fig.  1  and  similar  figures  have  been  used,  whenever  appro¬ 
priate,  to  single  out  those  distribution  functions  which  probably  will 
fit  well  to  the  observed  data.  This  is  done  by  examining  whether  the 
point  (bpb^)  plotted  in  the  P^ane  ls  inside  the  region  (or  close 

to  the  point  or  the  curve)  associated  with  the  distribution  function  for 
which  the  goodness-of-f it  is  to  be  considered, 

2.3  Test  of  Goodness-of-Fit 

Although  there  are  a  number  of  possible  ways  in  which  a  test  can  be 
performed  on  goodness-of-f it  (for  example,  the  y2  test  and  the  Kolmogorov- 
Smirnov  test),  the  w2  method  [3]  is  chosen  for  the  present  investigation. 
Let  F(x)  =  distribution  to  be  tested  and  Fn (x )  =  empirical  distribution 

based  on  a  sample  of  size  n,  and  then  form  a  statistic 

00 

nw2  =  n  /  [F  (x)  -  F(x)]2dF(x>  (5) 

—  CO 

This  statistic  is  distribution-free  and  some  of  the  percentiles  of  its 
asymptotic 'distribution  are  listed  below  [3]  and  are  also  plotted  in  Fig. 
2. 


Table  2  The  w2  Test 


P(noj2  ^  a} 

0.80 

0.85 

0.90 

0.95 

0.98 

0.99 

a 

0.241 

0.285 

0.347 

0.461 

0.620 

0.743 

For  practical  computations,  the  following  simpler  form  of  nco2  has  been  used. 
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Fig.  1  Regions  of  various  distributions  in  the  (g  , 0  )  plane 


Cumulative  Probability  P{nw 


Realization  of  nw2 

n 


Fig.  2  Probability  Distribution  of  Statistic  nw2 

n 
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nu>2  =  l/(12n)  +  l  [(21  -  l)/(2n)  -  F(x.)]2 
n  1=1  1 

=  l/(12n)  +  l  [(21  -  l)/(2n)  F(xm)]2  (6) 

1=1 

where  x^j  (1=1 ,2, . . .  ,n)  are  the  observations  arranged  In  ascending  order; 

Is  the  smallest,  x^j  the  second  smallest,  ...,  x^  the  largest. 

The  preceding  table  indicates  that  the  proposed  distribution  should 
be  rejected  if  the  corresponding  nw2  value  is  larger  than  0.461  (0.347) 
for  a  significance  level  of  5%  (10%). 

2.4  Statistics  of  Observed  Data 

There  are  three  sets  of  observed  data:  (i)  XQPF,  ( 1 i )  XWPF  and  (iii) 

WPF,  where  X  specif ies  "Load  Transfer",  W  specifies  "Winslow  Drilled",  P  speci¬ 
fies  "Proper  Technique,"  F  specifies  "Fighter  Spectrum"  and  Q  specifies 

-3 

"Quackenbush  Drill  and  Ream."  All  data  values  are  given  in  mils  (10 
inch) . 

The  observed  data,  XQPF,  XWPF  and  WPF  are  arranged  in  ascending  order 
as  shown  in  Table  1  with  the  cumulative  probability  or  the  plotting  posi¬ 
tion  defined  by 

F(x(.})  =  1/(n  +  1)  (7) 

Estimates  of  some  fundamental  statistics  are  also  listed  in  Table  1.  Sym¬ 
bols  used  therein  signify  the  following:  MEAN  =  mj,  AM2  =  m^,  AM3  =  m^, 

AM 4  =  m^,  STDV  =  s  (unbiased  standard  deviation),  SQRT(Bl)  =  v^,  B2  =  b2 
and  B1  =  b^. 
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SECTION  III 


FITTING  TO  THE  JOHNSON 
DISTRIBUTION  FAMILY 

3.1  Method  of  Translation 

The  Johnson  distribution  family  was  derived  by  Johnson  [2]  with  the 
aid  of  a  method  of  translation  which  takes  advantage  of  possible  transfor¬ 
mations  of  non-Gausslan  random  variables  into  Gaussian  (normal)  variables. 
The  method  is  outlined  below. 

We  say  that  a  random  variable  X  has  been  transformed  to  the  normality 
if  a  function  G(-)  transforms  X  into  the  standardized  normal  variable  Z. 


Z  =  G ( X )  (8) 

with  the  density  function  of  Z  being 

f z ( z )  =  <f)(z)  =  1/^ir  exp(-z2/2)  (9) 

Such  a  transformation  can  be  performed  in  two  steps.  First,  we  perform  a 
linear  transformation  of  X  into  Y  such  that 

Y  =  (X  -  e)/A  (10) 

with  A  being  positive  and  then  transform  Y  into  the  standardized  normal 
variable  Z  by 

Z  =  G(X)  =  y  +  ng(Y)  (11) 


where  n  is  assumed  to  be  positive.  The  density  function  of  Y  can  then  be 
shown  to  be 

fy(y)  =  ng'(y)fz(y  +  ng(y)) 

=  ng'(y)/v#nr  exp[-{y  +  ng (y) }2/2]  (12) 

where  g(y)  has  been  assumed  to  be  a  non-decreasing  function  of  y  and  g'(y) 
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=  dg(y)/dy.  We  note  that  there  are  four  parameters  e,  A,  y  and  n  involved 
in  the  transformation  Z  =  G(X). 

Since  X  and  Y  are  linearly  related  through  Eq.  10,  it  is  easy  to  show 
that  the  density  function  of  X  is  given  by 

fx(x)  =  n/(v&)g'(^-)  exp[-{y  +  ng(^-)>2/2]  (13) 

The  transformation  Z  =  G(X)  =  y  +  ng(Y)  and  the  relationship  between  f^(z) 
and  f  (x)  are  schematically  illustrated  in  Fig.  3.  Beyond  the  obvious  fact 
that  the  analytical  form  of  g ( * )  precisely  determines  the  distribution  func¬ 
tion  of  Y  when  n  =  1  and  y  =  0,  we  observe  from  Fig.  3  that  the  skewness 
and  the  kurtosis  of  the  distribution  of  Y  are  greatly  influenced  by  y  and 
n  respectively. 


3.2  The  Johnson  Distribution  Family 

Different  distribution  functions  can  be  generated  by  using  different 
functions  for  g(y).  In  fact,  Johnson  [2]  proposed  three  different  types 
of  distribution  functions  referred  to  as  the  Johnson  SL>  Sg  and  Sy  distri¬ 
butions  by  respectively  employing  the  following  three  functions  for  g (y ) : 

(a)  g(y)  =  Jtn(y)  for  y  >  0  (14) 

(b)  g(y)  =  £n(y^)  =  2tinh’:1  (2y  -  1)  for  0  <  y  <  1  (15) 

(c)  g(y)  =  sin  h_1(y)  =  £n{y  +  /y2 '+  1}  for  -°°  <  y  <  °°  (16) 

3.2.1  The  Johnson  SL  Family 

The  density  function  of  the  Johnson  family  can  be  defined  with  the 
aid  of  Eqs.  13  and  14  as 

fS|_(x)  =  ti/[Y^tt(x  -  e)]  exp{-[y  +  n£n(^-)]z/2}  for  x  >  e 

(17) 

where 


n  >  0,  jy|  <  °°,  X  >  0,  |e 
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<  OO 


y  is  defined  by 
0  =  y  +  nf(y) 


+ 

II 


c  to 
o  c 

»r—  fO 

•*->  >- 
o  I— 


=1  aj-— - 
u_  x:  >- 

4->  r-< 
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Fig.  3  Schematic  Representation  of  the  Method  of  Translation 


Setting 


y*  =  y  -  n£nA  (18) 

we  can  rewrite  Eq.  17  In  the  following  form: 

fSL(x)  =  n / [v^tt  (x  -  e)]  exp{ -%n2 [y*/n  +  £n(x  -  e)]2}, 

x  >  e  (19) 

where 


n  >  0,  | y* !  <  00 .  ! e |  <  00 

We  can  recognize  Eq.  19  as  the  three-parameter  log-normal  distribution.  In¬ 
deed,  introducing  a  and  y  so  that 

n  =  1/a  and  y*  =  -y/a  (20) 

we  can  derive  the  familiar  form  of  the  three-parameter  log-normal  density 
function : 


fSL(x)  =  l/[v^7ra  (x  -  e)]  exp{-Un(x  -  e)  -  y]2/(2a2)}  (21) 


With  the  form  of  the  transformation  function  given  in  Eq.  14,  we  can 
show  that  the  i-th  moment  about  the  origin  yl  of  variable  Y  is  given  as 

00 

yj(y)  =  l/( v^tt)  /  exp{i (z  -  y)/n>  exp(-z2/2)dz 
-00 

=  exp{(i/n)2/2  -  y(i/n)}  (22) 

It  then  follows  that  8^  (square  of  skewness)  and  ^  (peakedness)  are  given 
by 


8j  =  (w  -  l)(w  +  2) 2 

8^  =  3  +  (w  -  l)(w3  +  3w2  +  6w  +  6)  ^ 


(23) 


— Cr 

where  w  =  e  with  a  indicating  the  standard  deviation  of  the  log-normal 
distribution.  When  this  relationship  is  plotted  on  the  8j ~8g  plane,  we 
obtain  a  curve  indicating  those  values  of  8j  and  82  that  represent  the 
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log-normal  distribution  as  shown  In  Fig.  4. 

3.2.2  The  Johnson  Sg  Family 

On  the  basis  of  the  function  g(*)  defined  by  Eq.  15,  we  can  construct 
f$B(x)  =  1  Iftn  •  A/{(x  -  e)(A  -  x  +  e)}  exp{-[y  +  n^n(^-~|£)] 2/2} 

for  e  <  x  <  e  +  \  (24) 

where 


n  >  0,  |y|  <  oo,  \  >  0,  jej  <  00 

The  probability  distributions  with  the  density  function  given  by  Eq.  24  are 
said  to  belong  to  the  Johnson  Sg  family.  These  distributions  involve  four 
independent  parameters  and  consequently  Bj  and  B2  that  represent  this  fam¬ 
ily  of  distributions  can  take  on  those  values  within  the  domain  designated 
by  "Johnson  Sg  distribution"  in  Fig.  4.  Indeed,  this  domain  is  bound  by 
the  curve  representing  the  log-normal  distribution  and  a  straight  line  B2- 
Bj  -  1  =  0.  Above  this  straight  line  is  the  domain  representing  those  val¬ 
ues  of  Bj  and  b2  that  are  not  realizable. 

It  follows  from  Eq.  12  that  the  density  function  of  Y  =  (X  -  e)/\  is 
given  by 

fYSB^  =  n/[v^y(l  -  y)l  exp{-[y  +  n^n(y^)]  2/2} 

for  0  <  y  <  1 

It  also  follows  from  Eq.  11  that  Y  is  expressed  in  terms  of  Z  as 

Y  =  {1  +  exp [- ( Z  -  yj/q]}'1 

Hence,  the  median  value  y  of  Y  is 

y  =  [1  +  expC+y/n)]”1  (27) 

The  necessary  and  sufficient  conditions  for  bimodality  irrespective  of  the 
sign  of  y  are 


(25) 

(26) 
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Standardized  Measure  of  Peakedness 


C\J 

c a 


Measure  of  Skewness 


Fig.  4  Estimated  Result  of  f°r  Data  xQpp>  XWPF 

and  WPF  (Johnson  Distribution  Family) 
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n  <  l/*/2,  | y j  <  {l/r\)A  -  2n2  -  2ntanh”^  A  -  2nz 


(28) 


3.2.3  The  Johnson  S[J  Family 

Use  of  Eq.  16  in  Eq.  12  immediately  results  in  the  density  function 

of  Y: 

fYSlJ(y)  =  n/i^ir(y2  +  1)  exp{-[y  +  r)S,n(y  +/y2  +  l)]2/2} 

for  -00  <  y  <  co  (29) 

Similarly,  use  of  Eq.  16  in  Eq.  13  produces  the  density  function  of  X: 

fsu(x)  =  n /v4t{(x  -  e)2  +  X2} 

x  exp{-[y  +  n£n{(x  -  e)/X  +  /[x.  -  e)?/X2  +  1}] 2} 

for  -oo  <  x  <  °°  (30) 


where 


n  >  0,  -oo<y<oo)  x  >0,  -oo  <  e  <  oo 

The  probability  distributions  having  the  density  function  given  by  Eq.  30 
are  said  to  belong  to  the  Johnson  family.  As  in  the  case  of  the  John¬ 
son  Sg  family,  the  values  which  and  82  of  the  distributions  of  this 
family  can  take  are  confined  in  a  domain  below  the  log-normal  curve  as 
designated  by  "Johnson  Su  distribution"  in  the  gj-62  plane  in  Fig.  4. 

3.3  Fitting  Data  to  Distributions  of  the  Johnson  Family 

Observed  experimental  data  can  be  fitted  to  the  distribution  functions 
of  the'  Johnson  family  by  proceeding  with  the  following  steps: 

(1)  Determine  which  of  the  three  distribution  families  is  to  be  used. 

(2)  Estimate  the  parameters  of  the  selected  distribution. 

(3)  Compute  the  expected  cumulative  frequencies  of  the  fitted  distri¬ 
bution. 

(4)  Perform  a  goodness-of-f it  test  using  the  co2  method. 

The  objective  of  Step  (1)  can  be  accomplished  by  plotting  on  the  3j-82 
plane  the  estimates  bj  and  b2  of  and  82, respectively evaluated  with  the 
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aid  of  Eq.  4.  Observe  whether  or  not  the  point  Is  in  the  SD  domain, 

b 

Sjj  domain,  or  close  to  or  on  the  Sl  curve.  It  appears,  prudent  to 
assume  that  the  data  fit  the  Sg  (S^)  distributions  if  the  point  ( b-^ »b2 ) 

falls  in  the  Sg  (S^)  domain,  and  that  the  log-normal  distribution  is  a  like¬ 
ly  candidate  if  the  point  falls  on  or  close  to  the  SL  curve.  Indeed,  for 
the  observed  data  listed  In  Table  1,  all  three  points  representing  (bj.b,,) 
for  XQPF,  XWPF  and  WPF  fall  In  the  Sg  domain.  Therefore,  these  data  are 
expected  to  fit  well  to  the  Sg  distribution.  Also,  while  the  three  points 
are  not  particularly  close  to  the  SL  curve,  their  general  proximity 
to  the  curve  suggests  that  reasonable  fits  may  be  expected  between  the  da¬ 
ta  and  the  three-parameter  log-normal  distribution.  In  fact,  the  data  plot¬ 
ted  on  the  log-normal  probability  paper  (assuming  that  the  minimum  flaw 
size  =  0)  in  Fig.  5  suggests  that  a  more  than  satisfactory  fit  may  be  ob¬ 
served  particularly  for  XWPF  and  WPF  If  the  three-parameter  log-normal  dis¬ 
tribution  is  used. 

For  the  reasons  described  above  and  for  a  greater  familiarity  with  the 
log-normal  distribution  on  the  part  of  engineers,  the  Johnson  SL  (three- 
parameter  log-normal)  distribution  is  considered  first  for  the  purpose  of 
fitting  the  observed  data  in  Table  1.  Then,  a  fitting  procedure  will  be 
described  for  the  Johnson  Sg  distribution. 

3.3.1  The  Johnson  (Three-Parameter  Log-Normal)  Distribution 

Expressing  the  log-normal  probability  density  function  in  the  form  of 
Eq.  21  with  parameters  y  ,  o  and  e  ,  we  can  establish  the  parameter  estimation 
procedure  in  one  of  the  following  manners,  depending  cn  whether  or  not  the 
location  parameter  e  is  kncwn.  The  number  of  unknown  parameters  is  equal 
to  two  when  e  is  known;  otherwise  it  is  equal  to  three. 

(a)  When  e  is  assumed  to  be  known 

Since  )ln(X  -  e)  is  a  normally  distributed  random  variable  with  mean  y 

A  A 

and  standard  deviation  a,  the  estimates  y  and  a  of  y  anda  ,  respectively,  can 
be  obtained  in  a  manner  analogous  to  that  for  estimating  the  parameters  of 
the  normal  distribution: 
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Cumulative  Probability  (%) 


Equivalent  Initial  Flaw  Size  x  (mils) 


Fig.  5  Log-Normal  Probability  Plots  of  Data  Sets  XQPF,  XWPF  and  WPF 
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(31) 


y  =  (1/n)  l  x* 

1=1  1 

o  =  [(  l  (x*  -  y)2}/^5*  =  t{  1  x*2}/n  -  (y)2]'1 
1=1  1  1=1  1 


(32) 


or 

s  =  [(  l  (xt  -  y)2}/(n  -  l)]'2  =  /n/(n  -  1)  c  (33) 

1=1  1 

where  x|  =  5,n(xi-  -  e)  and  s  is  the  unbiased  estimate  of  standard  deviation. 

The  estimates  of  n  and  y*  in  Eq.  20  are  then  obtained  as 

A  A  A 

rij  =  1/a  or  n2  =  1/s  (34) 

and 

A  A  A  A  A 

Yj  =  -v/a  or  y*  =  -y/s  (35) 

Because  of  the  general  theoretical  advantage,  s  in  Eq.  33  is  used  more  frequent¬ 
ly  than  a  in  Eq.  32  for  an  estimate  of  the  standard  deviation.  For  a  large 
value  of  n,  however,  there  is  little  difference  between  these  two  as  Eq.  33 

A 

indicates.  In  the  present  section,  however,  both  a  and  s  are  used  for  com¬ 
parison  . 


A 

As  is  well  known,  the  estimate  y  in  Eq.  31,  when  realizations  x^.’s  are 
replaced  by  corresponding  random  variables  ' s ,  becomes  an  estimator  and 
will  vary  from  sample  to  sample.  Therefore,  it  is  standard  practice  that 

A 

the  adequacy  of  y  as  an  estimate  for  the  unknown  parameter  y  is  indicated 
in  terms  of  the  confidence  interval  given  by 


yg 


=  y  ±  (t 


l-a,n-l 


)s/>/n 


(36) 


A 

where  y  and  s  are  obtained  from  Eqs.  31  and  33,  respectively,  and  t,  , 

X jP ”1 

is  the  two-sided  100(l-a)  percentile  of  the  Student's  t  distribution  of  (n-1) 
degrees -of -freedom.  The  quantity  (1-a)  is  referred  to  as  the  confidence  le¬ 
vel  and  has  the  following  significance. 
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)  s/i/n}  =  1  -  a 


(37) 


^  -  (tl-a,n-l)s/^  ^  ^  +  (tl-c.,n-l 

From  the  (two-sided)  100(l-a)%  confidence  interval  of  y  =  E[£n(X  -  e)] 
given  in  Eq.  36,  we  can  obtain  the  lower  bound  XL  and  upper  bound  Xy  of  the 
corresponding  interval  for  X  from 

yL  =  i,n(XL  -  e),  yy  =  £n(Xy  -  e)  (38) 

as 


XL  =  e  +  exp{y  -  (tj  )sMf} 

(39) 

A 

xy  =  e  +  exp{y/+  (ti_afn-l^s/^ 

(40) 

(b)  When  e  is  assumed  to  be  unknown 


It  follows  from  Eqs.  11,  14,  and  18  that 


Z  =  y*  +  n£n(X  -  e)  (41) 


where  Z  is  the  standardized  normal  variable.  The  three  parameters  y*,  n 
and  e  in  Eq.  41  can  be  estimated  in  the  following  manner  as  suggested  by  Hahn 
and  Shapiro  [4]:  choose  three  probability  values  p,  q  and  r,  find  A  =  $(zA) 
where  $(•)  =  standardized  Gaussian  distribution  function,  with  A  =  p,  q  and 
t ,  estimate  such  that  P{X  x^}  =  A,  construct  three  equations  of  the 
form 


zA  =  Y*  +  n&n(xA  -  e)  (A  =  p,  q  and  r)  (42) 

A  A  A 

and  finally  solve  Eq.  42  for  the  estimates  y*,  n  and  e.  In  practice,  the 
quantity  x^  defined  above  is  estimated  with  the  aid  of  the  ordered  sample 

X(l)  <  x^)  <  - <  x^nj  where  x^  is  the  i-th  smallest  in  the  sample  of 

size  n  with  the  plotting  position  i/(l+n).  Indeed,  if  it  so  happens  that 


A  =  i/(l+n) 


(43) 


then 


x 


A 


(44) 
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while  if 


1/(1  +  n)  <  A  <  (1  +  1)/(1  +  n)  (45) 

then  by  interpolation 

XA  =  x(i )  +  {x(i+l )  "  x(1)^^  +  n^A  “  (46) 

In  the  present  study,  we  choose  p  =  1-a,  q  =  0.5  and  r  =  a  with  a  being 
i/(l+n)  (1=1 ,2, . . . ,8) .  Then,  it  is  a  relatively  simple  matter  to  derive 

A  A  A 

the  following  expressions  for  y*.  n  and  e: 


Y*  -  n£n{(l  -  e‘z  ^n)/(*Q  5  -  xj) 

(47) 

1  ■  z'/Jln{(x1.a  -  V5)/tx0.5  ' 

A  A  A 

(48) 

e  =  xQ  5  -  exp(-y*/n) 

(49) 

Obviously,  these  values  are  different  for  different  values  of  a.  In  the 

A 

analysis  that  follows  later,  we  use  as  our  estimates  those  values  of  y*, 

/N  /S 

n  and  e  that  produce  the  smallest  value  of  nw*.  In  Eqs.  47  and  48,  z'  = 
Zi  -  -z  represents  the  100  x  (l-a)-th  percentile  of  the  standardized 
Gaussian  distribution. 


(c)  Results  of  estimation 


First,  we  consider  the  case  where  the  location  parameter  e  is  assumed 
to  be  known.  In  this  case,  we  perform  the  estimation  presuming  that  e  is 


a  fraction  of  the  smallest  observation  x. 


(1 ) ’  £  =  X(l) 

9).  Using  Eqs.  31-35,  we  then  estimate  n  and  y*  based  on  both  o  and  s  for 


i/10  (1=0,1, 2, . . 


each  of  these  ten  different  values  of  e;  e  -  0,  O.lx 


(1)’ 


0.9x 


(1 )  * 


For  each  set  of  n  and  y*  thus  estimated,  Eq.  6  is  used  to  evaluate  nor, 

a  a  n 


and  we  choose  as  our  best  estimate  the  set  of  e,  n  and  y*  that  produce 
the  smallest  value  of  nw*.  Table  3  lists  estimated  parameters  (writ¬ 
ten  as  ETA  1),  y*  (GAMMA  1),  ^  (ETA  2),  y*  (GAMMA  2)  and  noo*  (NWN2)  for 
the  data  XQPF,  XWPF  and  WPF  under  ten  different  cases  of  e;  Case  1  for 
e  =  0,  Case  2  for  e  =  O.lx^y  ...,  Case  10  for  e  =  0.9x^.  The  nto* 
values  are  considerably  larger  for  XQPF  than  for  XWPF  and  WPF.  The  noo* 


20 


Table  3  Values  of  n.  Y  and  na£  as  Functions 

f  II  . 

of  e  (The  Johnson  SL  Distribution) 


FOR  DATA  SERIES  XQPF 


CASE 

E 

ETA  1 

GAMMA  1 

NWN2 

ETA2 

GAMMA2 

NwN? 

1 

.0000 

.736 

1.119 

1.630 

.726 

1.109 

1  .696 

2 

.0026 

.725 

1 .120 

1.673 

.715 

1.105 

1  .690 

3 

.0052 

.719 

1.120 

1.720 

.709 

1.105 

1.738 

9 

.0078 

.702 

1.120 

1.771 

.692 

1.105 

1.790 

5 

.0)09 

.689 

1.120 

1.828 

.680 

1.109 

1  .899 

6 

.0130 

.676 

1.118 

1.892 

.666 

1.103 

1.919 

7 

.0156 

.661 

1.115 

1.965 

.652 

1.100 

1.988 

6 

.0182 

.699 

1.110 

2.050 

.635 

1.095 

2.075 

9 

.0208 

.623 

1.101 

2.155 

.615 

1.086 

2.182 

10 

.0239 

.595 

1.083 

2.300 

.587 

1.069 

2.331 

FOR  DATA 

i  SERIES 

XrtPF 

CASE 

E 

ETA  1 

gammai 

NftN2 

ETA2 

GAMMA2 

NtoN2 

1 

.0000 

1  .671 

1  .905 

.951 

1.698 

1.879 

.955 

2 

.0093 

1.621 

1.909 

.971 

1.599 

1.879 

.975 

■  3 

.0186 

1.569 

1  .901 

.993 

1.598 

1.876 

.997 

a 

.0279 

1.516 

1.895 

.519 

1.995 

1.870 

.523 

5 

.0372 

1.960 

1.886 

.598 

1.990 

1.860 

.553 

6 

.0965 

1.902 

1.672 

.582 

1.383 

1.896 

.587 

7 

.0558 

1.390 

1.851 

.622 

1.321 

1.826 

.628 

8 

.0651 

1.271 

1.822 

.673 

1.259 

1.797 

.680 

9 

.0799 

1.191 

1.776 

.791 

1.175 

1.751 

.798 

10 

.0837 

1.085 

1.691 

.896 

1.070 

1.668 

.859 

FOR  DATA 

SERIES 

wPF 

CASE 

E 

ET  A  1 

GAMMA  1 

NWN2 

ETA2 

GAMMA2 

NWN2 

1 

.0000 

1.586 

.719 

.330 

1.565 

.709 

.333 

2 

.0190 

1.598 

.738 

.393 

1.527 

.728 

.397 

3 

.0280 

1.507 

.761 

.358 

1.987 

.751 

.362 

H 

.0920 

1.965 

.783 

.379 

1.995 

.772 

.378 

5 

.0560 

1  .920 

.802 

.393 

1.901 

.792 

.398 

6 

.0700 

1.372 

.820 

.915 

1.359 

.809 

.921 

7 

.0890 

1.320 

.639 

.993 

1.302 

.823 

.998 

8 

.0980 

1  .260 

.699 

.977 

1.293 

.833 

.983 

9 

.1120 

1.188 

.896 

.525 

1.172 

.835 

.532 

10 

.1260 

1.085 

.828 

,606 

1.071 

.817 

.619 
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values  for  the  last  two  sets  are  not  significantly  different.  Hence,  the 
three-parameter  log-normal  distribution  can  be  used  for  the  XWPF  data  with 
approximately  the  same  level  of  goodness-of-flt  as  for  the  WPF  data,  while  at 
a  considerably  less  satisfactory  level  for  the  XQFF  data,  a  trend  that  can 
easily  be  observed  fran  Fig.  5.  Table  3  also  shows  that  within  each  set  of 
data,  the  nut:  values  decrease  as  e  decreases  thus  indicating  the  choice  of 
e  =  0  together  with  the  corresponding  values  of  n  and  y*  as  statistically 
the  best.  The  value  of  nw2  associated  with  the  significance  level  of  5% 
is  0.461  from  Table  2.  Therefore,  observing  from  Table  3  that  the  nw2 
values  for  e  =  0  are  smaller  than  0.461  for  XWPF  and  WPF,  we  may  accept 
the  hypothesis  (with  a  significance  level  of  5%)  that  XWPF  and  WPF  data 
have  been  taken  from  the  three-parameter  log-normal  populations  with  e  = 

A  A  A  A 

0  and  corresponding  values  of  rij  and  y£  (or  and  y|).  However,  we  must 
reject  (with  the  same  level  of  significance)  the  hypothesis  that  XQPF  da¬ 
ta  have  been  taken  from  the  three-parameter  log-normal  populations  since 
the  smallest  nw2  value  associated  with  e  =  0  is  in  this  case  larger  than 
0.461.  It  is  of  Interest  to  note  that  if  the  significance  level  is  raised 
to  10%,  only  WPF  data  will  survive  the  test.  These  results  are  summarized 
In  Table  4. 


Figures  6-8  show  the  values  of  y£,  y|,  rij  and  rv>  as  functions  of  the 
location  parameter  e  for  XQPF,  XWPF  and  WPF  data,  respectively.  We  ob¬ 
serve  from  these  figures  that  the  difference  between  the  estimates  based 

A 

on  o  and  s  are  neglibibly  small.  Fig.  9  illustrates  how  the  values  of 

A  A  A 

ni  and  y|  compare  "data  set  by  data  set"  as  functions  of  e  when  a  is  used, 
while  Fig.  10  illustrates  the  same  when  s  is  used.  These  two  figures  al- 

A  A 

so  show  that  the  values  of  both  n  and  y*  are  generally  largest  for  XWPF, 

larger  for  WPF  and  smallest  for  XQPF  over  practically  the  entire  range 

of  e  considered.  Table  4  also  lists  the  upper  and  lower  bounds,  Xy  and 

X^,  corresponding  to  the  confidence  bounds  py  and  p^  with  a  confidence 

level  of  0.9,  while  Table  5  lists  the  bounds  corresponding  to  confidence 

levels  of  0.9,  0.95  and  0.99.  Fig.  11  plots  the  n«2  values  as  functions 

of  e  respectively  for  XQPF,  XWPF  and  WPF  and  in  essence  reiterates  the 

result  of  the  test  of  hypothesis  mentioned  earlier.  Note  that  the  nw2 
a  n 

values  based  on  a  and  s  are  indistinguishable  in  Fig.  11  for  the  same 
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Table  4  Best  Estimates  (The  Johnson  Distribution  with  e  =  0) 


Data 

e 

A 

n  . 

IH 

■9 

■91 

XQPF 

0 

.726 

1.10 

1.65 

.163 

.294 

XWPF 

0 

1.65 

1.88 

.455 

.281 

.364 

WPF 

0 

"1.57 

.704 

.333 

.557 

.730 

*  5%  (10%)  significance  level  =  .461  (.347) 

**  90%  confidence  bounds  In  mils 


Table  5  Lower  and  Upper  Bounds  XL  and  Xjj  In  mils  Corresponding  to 
,  and  py  for  Several  Confidence  Levels  1-a  (When  e  =  0) 


FOR  DATA  XQPF 

( N®37 ) 

1-ALPHA 

XL 

XU 

0.90 

.1627 

.2938 

0.95 

.1492 

•  3204 

0.99 

.1260 

.3793 

FOR  DATA  XWPF 

(N=37) 

1-ALPHA 

XL 

XU 

0.90 

.2808 

.3643 

0.95 

.2703 

.3785 

0.99 

.2509 

,4077 

t  FOR  DATA  WPF 

(N=38) 

1-ALPHA 

XL 

XU 

0.90 

.5571 

.7299 

0.95 

.5354 

.7595 

0.99 

.4957 

,8203 
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Variation  of  Estimated  Parameters  n  and  y*  for  Data  XQPF 
(When  e  is  Known) 
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Fig.  8  Variation  of  Estimated  Parameters  th  and  y*  for 'Data  WPF 

(When  e  is  Known) 
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Unbiased  Statistic  s  (When  e  is  Known) 


jCou  san[EA 


Fig.  11  Effect  of  the  Assumed  Location  Parameter  e  on  the 
Gpodness-of-Flt  Test  Statistic  nu>^ 


data  sets. 


We  now  turn  to  the  case  where  e  Is  assumed  to  be  unknown.  In  this  case, 

A  A  A 

Eqs.  47-49  are  used  to  estimate  e,  n  and  y*.  As  mentioned  earlier,  differ- 

AAA 

ent  sets  of  e,  n  and  y*  will  be  obtained  depending  on  the  values  of  a  to  be 
used  in  Eqs.  47  and  48.  In  the  present  study,  we  choose  the  following  eight 
different  values  for  a;  a  =  i/(l+n)  (i=l ,2 , . . .  ,8) .  Making  use  of  either 
Eqs.  43  and  44  of  Eqs.  45  and  46,  the  values  of  x  ,  x1  and  xQ  5  are  eval¬ 
uated  and  listed  in  Table  6  for  each  case  of  these  a  values;  Case  i  for  a 

A 

=  1/(1 +n ) .  Table  7  lists,  for  each  data  set,  the  estimated  parameters  n 

A  A 

as  ETA,  y*  as  GAM,  e  as  E  and  it  further  lists  the  value  of  nw2  computed 

A  A  A  n 

with  the  aid  of  Eq.  6.  The  values  of  n,  y*  and  e  for  XQPF,  XWPF  and  WPF 
are  plotted  respectively  in  Figs.  12-14  as  functions  of  the  probability 
level  a.  Fig.  15  plots  the  value  of  nw2  for  each  data  set  as  a  function  of 

a  and  shows  that,  as  in  the  case  of  Fig.  11,  the  values  are  largest  for 

XQPF,  larger  for  XWPF  and  smallest  for  WPF,  again  reflecting  the  degrees  of 
"goodness-of-fit"  observed  in  Fig.  5.  Fig.  15  further  shows  that  the  nw2 
values  assume  minimum  at  a  =  3/38  =  0.0789  for  XQPF,  at  a  =  1/38  =  0.0263 

A  A  A 

for  XWPF  and  at  a  =  5/39  =  0.1316.  The  set  of  e,  n  and  y*  corresponding 
to  each  of  these  a  values  is  chosen  as  the  best  estimate  for  the  respec¬ 
tive  data  set.  If  we  use  the  best  estimate  of  e  thus  obtained  in  place  of 

e  in  Eqs.  39  and  40,  the  lower  and  upper  bounds,  X^  and  X^,  will  result  as 
listed  in  Table  8.  In  evaluating  X^  and  X^  from  Eqs.  39  and  40,  y  and  s 

are  needed  and  they  are  computed  with  the  aid  of  Eqs.  31  and  33. 

It  is  concluded  from  Table  7  and  Fig.  15  that  with  the  significance 
level  of  5%  we  may  accept  the  hypothesis  that  the  WPF  data  have  been  taken 
from  the  three-parameter  log-normal  population  but  we  must  reject  the  oth¬ 
er  two  data  sets  as  taken  from  the  three-parameter  log-normal  populations. 

A  A  A 

For  the  WPF  data,  the  best  estimates  are  n  =  1.855,  y*  =  0.624  and  e  = 

-0.067.  Incidentally,  the  WPF  data  will  also  pass  the  hypothesis  testing 
under  the  significance  level  of  10%.  These  results  are  summarized  in 
Table  9. 
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Table  6  Values  of  x  5Q,  xfl  and  x1-0 
(The  Johnson  SL  Distribution) 


XGPF 

XI4PF 

HPF 

XC 

.50) 

.150 

.295 

.647 

CASE 

1 

XC 

.03) 

.026 

.093 

.140 

X( 

.97) 

7.700 

1.260 

3.830 

CASE 

2 

xc 

.05) 

.031 

.145 

.236 

X( 

.95) 

3.000 

.940 

2.730 

CASE 

3 

X( 

.08) 

.0«7 

.145 

.280 

xc 

.92) 

1.240 

.870 

1  .640 

CASE 

4 

xc 

.11) 

.050 

.160 

.280 

x  c 

.90) 

1.140 

.610 

1.490 

CASE 

5 

xc 

.13) 

.058 

.175 

.320 

xc 

.87) 

1.100 

.810 

1  .250 

CASE 

6 

xc 

.16) 

.059 

.180 

.367 

XC 

.85) 

1.090 

.700 

1.250 

CASE 

7 

XC 

.18) 

.060 

.180 

.367 

XC 

.82) 

1.090 

.612 

1.140 

CASE 

8 

xc 

.21) 

.063 

.190 

.420 

xc 

.79) 

1.090 

.540 

1.040 
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Table  7  Values  of  n.  I.  e  and  nw* 

(The  Johnson~S^  Distribution) 


DATA 

XQPF 

XKPF 

WPF 

El  A 

.<*72 

1.223 

1.061 

CASE  1 

GAM 

,977 

1.676 

.537 

E 

.024 

.041 

.044 

NwN2 

3.098 

.719 

.465 

ETA 

.510 

1.111 

1.006 

CASE  2 

GAM 

1.064 

1.613 

.673 

6 

.026 

.100 

.135 

NaN2 

2.898 

.838 

.517 

El  A 

.599 

1.051 

1.433 

CASE  3 

GAM 

1.301 

1.676 

.775 

E 

.036 

.092 

.065 

NaN2 

2.408 

.872 

.340 

ETA 

.546 

.935 

1.524 

CASE  4 

GAM 

1  .200 

1  .586 

.656 

E 

.039 

.112 

-.003 

NWN2 

2.600 

.999 

.322 

ETA 

.479 

.768 

1.855 

CASE  5 

GAM 

1.095 

1.425 

.624 

E 

.048 

.139 

-.067 

NAN2 

2.717 

1.190 

.278 

El  A 

.430 

.797 

1.330 

CASE  6 

GAM 

.966 

1,457 

.863 

E 

.049 

.134 

.124 

NAN? 

2.940 

1.178 

.359 

ETA 

.383 

.887 

1.622 

CASE  7 

GAM 

'  .885 

1.519 

.703 

E 

.050 

.115 

-.001 

NWN2 

3.299 

1.059 

.306 

ETA 

.338 

.950 

1.500 

CASE  8 

GAM 

.793 

1.609 

.932 

E 

.054 

.111 

.110 

NwN2 

4.077 

.983 

.322 
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2.00 


Fig.  12  Variation  of  Estimated  Parameters  n»  Y*  and  e  for  Data  XQPF 

(When  e  is  Unknown) 
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Fig.  13  Variation  of  Estimated  Parameters  n,  Y*  and  e  for  Data  XWPF 

(When  e  i$  Unknown) 
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Variation  of  Estimated  Parameters  n,  T*  and  e  for  Data  WPFr 
(When  e  is  Unknown) 
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Fig.  15  Effect  of  the  Probability  Level  a  on  the  Goodness-of-Fit 

Test  Statistic 


Table  8  Lower  and  Upper  Bounds  XL  and  In  mils  Corresponding  to 
and  Uy  for  Several  Confidence  Levels  1-a  (The  Johnson 
Distribution;  When  e  Is  Assumed  to  be  Unknown) 


FOR  DATA  XQPF 

(N=37) 

1-ALPHA 

XL 

XU 

0.90 

.1655 

.2988 

0.95 

.1526 

.3276 

0.99 

.1314 

.3929 

FOR  DATA  XWPF 

(Ns37) 

1-ALPHA 

XL 

XU 

0.90 

.2737 

.3546 

0.95 

.2637 

.3687 

0.99 

.2455 

.3978 

FOR  DATA  WPF 

(N=38) 

1-ALPHA 

XL 

XU 

0.90 

.5673 

.7425 

0.95 

.5450 

.7721 

0.99 

.5039 

.8325 

Table  9  Best  Estimates  (The  Johnson  SL  Distribution  With  Unknown  e) 


Data 

A 

n 

A 

Y* 

A 

e 

n“S  * 

v  ** 

*L 

y  ** 

*U 

XQPF 

.599 

1.30 

.036 

2.41 

.166 

.299 

XWPF 

1.22 

1.68 

.041 

.719 

.274 

.355 

WPF 

1.86 

.624 

-.067 

.278 

.567 

.743 

5 %  (10%)  significance  level  =  .461  (.347) 
**  90%  confidence  bounds 


3.3.2  The  Johnson  $b  Distribution 

The  Johnson  SD  variable  Is  limited  between  the  lower  bound  e  and  the 

D 

upper  bound  e  +  A  as  Indicated  In  Eq.  24.  Therefore,  the  following  four 
possibilities  arise  with  respect  to  the  upper  and  lower  bounds:  (i)  Both 
bounds  are  known,  ( i 1 )  only  the  upper  bound  is  known,  (iii)  only  the 
lower  bound  is  known  and  (1v)  neither  bound  is  known.  In  the  present 
study,  we  formulate  the  procedures  of  parameter  estimation  assuming  that 
either  case  (iii)  or  case  (iv)  will  prevail,  although  actual  estimations 
are  performed  only  for  case  (iii).  It  is  noted  in  this  connection  that 
the  flaw  size  can  never  be  negative  and  therefore  the  lower  bound  may  be 
assumed  to  be  zero,  an  assumption  that  generally  produces  a  conservative 
result.  While  such  an  assumption  offers  a  considerable  mathematical  con¬ 
venience,  there  is  no  definite  reason,  physically  or  otherwise,  to  believe 
that  the  lower  bound  of  the  initial  flaw  size  distribution  must  be  equal 
to  zero. 


(a)  When  e  is  assumed  to  be  known 


In  this  case,  we  estimate  the  parameters  A,  n  and  y  assuming  that  e  = 

•  i/10  (1=0,1, . . .  ,9)  as  was  done  when  dealing  with  the  SL  distribution. 
Following  then  the  same  method  that  produced  Eqs.  47-49  from  Eq.  42,  we  ob- 

A  A  A 

tain  the  three  equations  below  for  the  estimates  A,  n  and  y. 


(x0.5  -  e)(e  +  X  '  xl-a' 1 


-Z  /Zi  I 
*1  a'  1  -a ' 


[(e  +  A  -  xQ  5)(xl  a,  -  e)J 


<*a  -  e)<E  +  X  '  x0.5) 
- - - 

(e  +  A  -  x  )(xQ  5  -  e) 


(50) 


A 


-  z  )/£n 

a 


(x,  ,  -  e)(e  +  A  -  x  ) 

1-a  a 

/\ 

(x  -  e)(e  +  A  -  x,  , ) 
'a  1-a 


(51) 


Y 


z,  ,  -  n  £n 
1-a 


x,  i  -  e 
1-a1 


e  +  A 


xi  i 

1-a 


(52) 


where,  as  before,  z^  and  x^  are  such  that  z^  =  $(A)  and  P{X  £  x^}  =  A  with 
$(•)  denoting  the  standardized  normal  distribution  function.  If  a  /  a' , 
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these  three  equations  must  be  solved  by  trial  and  error.  If,  however,  we 

take  an  Identical  value  for  a  and  a*  (a  =  a'),  then  Eq,  50  can  Immediately 

* 

be  solved  for  X  producing 

X  =  (*0.5  -  e>{(x0.5  -  e><xa  ■  e)  +  <x0.5  ‘  e>'xl-c,  '  e) 

-2(xa  -  €)(xl  o  -  e)>/{(x0  5  -  e)2  -  (xa  -  eJtXj^  -  e) 


A  A 

(53) 

The  estimates 

n  and  y  are  then  given 

by 

n  = 

2zl-a/tS,n^xl-ct  '  e^e 

A  A 

+  X  -  x  )}  -  £n{(x  +  e)(e  +  X  ■ 

a  a  ' ' 

'  xl-a 

/\ 

A 

(54) 

Y  = 

2l-a  -  ,'{*n<xl-a  -  • 

-  £n(e  +  X  -  x,  ) 

1-a 

(55) 

We  point  out  again  that  the  estimated  values  of  these  parameters  depend  on 
the  value  of  a  and  that  the  set  of  estimated  values  producing  the  smallest 
nco*  value  will  be  considered  as  the  best  estimate. 

(b)  When  neither  upper  nor  lower  bound  is  known 

It  follows  from  Eqs.  11,  15  and  18  that  the  Johnson  Sg  variable  X  and 
the  standardized  normal  variable  Z  are  related  by 

Z  =  y  +  n  £n{(X  -  e)/(e  +  X  -  X)}  (56) 

With  the  aid  of  z.  and  x.,  we  then  derive  the  following  four  equations 

n  Ma  a  a  a 

from  which  the  estimates  e,  X,  y  and  n  can  be  solved. 

A  /V  A  A  A 

zA  ■  Y  +  n  £n{(xA  -  e)/(e  +  X  -  xA)  (A=p,q,r  and  u)  (57) 

The  similarity  between  Eqs.  56  and  57  and  Eqs.  41  and  42  is  obvious.  As 
before,  depending  on  the  probability  levels  p,  q,  r  and  u  to  be  used,  dif¬ 
ferent  sets  of  estimates  will  be  obtained. 

However,  we  have  not  pursued  this  avenue  of  investigation  since  the 
preliminary  result  indicated  that  the  fit  of  the  observed  data  to  the  SD 

D 

distribution  with  four  unknown  parameters  would  probably  not  be  particu- 
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larly  better.  If  not  worse,  than  the  fit  to  the  same  distribution  with  three 
unknown  parameters.  In  this  connection,  we  recall  the  distribution  re¬ 
sult  where  considering  e  as  an  unknown  does  not  really  improve  the  goodness- 
of-fit. 


(c)  Confidence  Interval 

AAA  A 

Using  the  estimated  parameters  X,  y»  n  and  e  (e  in  case  it  is  assumed 
to  be  known),  write  the  following  equation 

Z  =  y  +  ri  £n{(X  -  e)/(e  +  X  -  X)}  (58) 


Since  Z  is  the  standardized  normal  variable,  the  unbiased  estimates  of  its 
expected  value  y  and  standard  deviation  o  are  respectively  given  by 


n 

V  =  l  z-/n  (59) 

i=l  1 


s  =  [{  l  (z. 
1=1  1 


-  ^)2}/(n  -  l)]h 


(60) 


where 

A  A  H  AAA 

z.  =  y  +  n  x  l  fcn{(x.  -  e)/(e  +  X  -  x . )}  (61) 

1  1=1  1  1 

/\ 

The  upper  bound  yy  and  the  lower  bound  y^  of  the  confidence  interval  of  y 
can  then  be  established  on  the  basis  of  the  Student's  t  distribution  as  in  the 
case  of  the  distribution.  Indeed,  they  are  also  given  by  Eq.  36.  The 
corresponding  bounds  X^  and  Xy  are  obtained  from 


XL  =  e  +  X  exp{(yL  -  y)/r\}/  [1  +  exp{(yL  -  Y)/n>l  (62) 


A  A  A  A  A  A 

Xy  =  e  +  X  exp{(yy  -  Y)/n)/[l  +  exp{(yy  -  y)/o)1 


(63) 


(d)  Results  of  estimation 

Dealing  with  the  case  where  e  is  assumed  to  be  known,  we  list  in  Table 
10  the  estimated  values  of  X,  n  and  y  of  the  Sg  distribution.  In  this  table. 
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Table  10  Results  of  Estimation  for  the  Johnson  Sg 
Distribution  When  e  Is  Assumed  to  be  Known 


FOR  DATA  SERIES  XWPF 


CASE 

E 

ET  A 1 

GAMMAi 

RAM 

NWN2 

1 

.0000 

•3^'^ 

T,  || 

_ _ *. 

2 

.0093 

-  . 

_ 

_ 

3 

.0186 

— 

4 

.0279 

-  ' 

_ 

.  .  ~ 

5 

.0372 

1.354 

4.774 

9.009 

.664 

6 

.0465 

1.260 

3.454 

4.103 

.669 

7 

.0558 

1.171 

2.703 

2.644 

.675 

8 

.0651 

1.088 

2.182 

1.939 

.683 

9 

.0744 

1.009 

1.791 

1.522 

.692 

10 

.0837 

.935 

1.484 

1  .244 

.703 

FOR  DATA 

SERIFS 

WPF 

CASE 

E 

ETA  1 

GAMMAI 

RAM 

NWN? 

1 

.0000 

1.512 

3.568 

7.503 

.283 

2 

.0140 

1.447 

3.143 

6.187 

.285 

3 

.0280 

1.385 

2.787 

5.253 

.287 

4 

.0420 

1.324 

2.484 

4.555 

.289 

5 

.0560 

1.265 

2.221 

4.012 

.291 

6 

.0700 

1.281 

2.832 

5.838 

.316 

7 

.0840 

1.221 

2.480 

4.852 

.319 

8 

.0980 

1.163 

2.185 

4.140 

.322 

9 

.1120 

1.486 

8.033 

119.569 

.322 

10 

.1260 

1.407 

4.907 

17.577 

.323 
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Case  1  indicates  the  use  of  e  =  xm  *  (i-l)/10,  and  E  is  written  for  e, 

ETA  1  for  n»  GAMMA  1  for  y,  RAM  for  A  and  NWN2  for  nw2 .  As  mentioned 

n 

earlier,  the  Sg  variable  Is  limited  between  e  and  e  +  X  (A  >  0) .  We 
have  found,  however,  that  the  estimation  procedure  used  here  produces 

A 

negative  values  for  the  estimate  A  depending  on  the  assumed  values  of  e. 

a 

For  example,  the  values  of  the  estimate  A  have  turned  out  to  be  nega¬ 
tive  for  all  ten  cases  of  e  with  respect  to  XQPF  data  and  for  the  first 
four  cases  with  respect  to  XWPF  data.  If  an  assumed  value  of  e,  say  en, 
produces  a  negative  A,  we  Interpret  that  to  be  an  indication  of  the  un- 
acceptability  of  the  Sg  distribution  with  e  =  £q  for  the  data  considered. 
Therefore,  the  results  of  the  estimation  for  XQPF  data  have  not  been 
listed  in  Table  10.  Also,  the  results  for  the  first  four  cases  of  e  val¬ 
ues  with  respect  to  XWPF  data  have  been  indicated  by  bars  (-). 

Based  on  the  nw2  values  listed  in  Table  10,  we  conclude  that  the  best 
n 

estimates  are  obtained  from  Case  5  for  XWPF: 

e  =  0.0372,  A  =  9.009,  ^  =  1.354,  y  =  4.774 
and  from  Case  1  for  WPF : 

e  =  0.000,  A  =  7.503,  r\  =  1.512,  y  =  3.568 

A  A  A 

Figures  16  and  17  plot  the  values  of  estimates  A,  n  and  y  as  functions 
of  e,  respectively  for  XWPF  and  WPF.  Fig.  18  plots  nio2  values  as  functions 
of  e  for  both  XWPF  and  WPF.  This  figure  shows  that  the  Sg  distribution  can 
be  accepted  for  WPF  data  but  not  for  XWPF  data  if  the  significance  level  of 
5 %  is  assumed.  Table  11  lists  the  values  of  X^  and  Xy  for  XWPF  and  WPF  with 
the  aid  of  Eqs.  62  and  63.  These  results  are  summarized  in  Table  12. 
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Johnson  Sn  Dist 
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Fig,  16  Variation  of  Estimated  Parameters  A,  n  and  y  in  the  Johnson 
SB  Distribution  for  Data  XWPF  (When  e  is  Known) 


10  00 


< 
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nation  of  Estimated  Parameters  X,  n  and  y  in  the  Johnson 
Sb  Distribution  for  Data  WPF  (When  e  is  Known) 


(The  Johnson  SR  Distribution) 
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Fig.  18  Effect  of  the  Assumed  Parameter  e  on  the -Goodness-of-Fit 
Test  Statistic  n^^  for  the  Johnson  S„  Distribution 


Table  11  Lower  and  Upper  Bounds  X^  and  In  mils  Corresponding  to 
and  for  Several  Confidence  Levels  1-a  (The  Johnson 
SR  Distribution;  When  e  Is  Assumed  to  be  Known) 


FOR  DATA  XQPF 

(N*37) 

1-ALPHA 

XL 

XU 

0.90 

.9567 

1.0314 

0.95 

.9485 

1.0925 

0.99 

.9289 

1,0645 

FOR  DATA  XWPF 

CN=37) 

1-ALPHA 

XL 

XU 

0.90 

.5270 

.8752 

0.95 

.4892 

.9429 

0.99 

.4232 

1.0898 

FOR  DATA  WPF 

(N=38) 

1-ALPHA 

XL 

XU 

0.90 

.5613 

.9586 

0.95 

.5189 

1.0370 

0.99 

.4455 

1.2079 

Table  12  Best  Estimates  (The  Johnson  Sg  Distribution) 


Data 

A 

X 

A 

0 

A 

Y 

e 

n<* 

Y  ** 

*L 

Y  ** 

XQPF 

A 

X  Is  negative;  unacceptable 

XWPF 

9.01 

1.35 

4.77 

.0372 

.664 

.527 

.875 

WPF 

7.50 

1.51 

3.57 

.00 

.283 

.561 

.959 

*  5%  (10%)  significance  level  =  .461  (.347) 

**  9025  confidence  bounds 
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SECTION  IV 


FITTING  TO  THE  WEIBULL 
DISTRIBUTION  FUNCTION 


An  attempt  is  made  to  fit  the  data  shown  in  Table  1  to  the  Weibull  dis¬ 
tribution.  For  the  present  investigation,  we  consider  the  following  three- 
parameter  Weibull  distribution. 

Fx(x)  =  1  -  exp{-[(x  -  eJ/X]11}  for  x  >  e  .  (64) 
where  n  =  shape  parameter,  X  =  scale  parameter  and  e  =  location  parameter. 


It  is  often  a  practice  to  use  the  maximum  likelihood  method  to  esti¬ 
mate  these  parameters.  Its  use  for  the  three-parameter  Weibull  distribu¬ 
tion,  however,  not  only  creates  numerical  problems  requiring  possibly  high¬ 
ly  expensive  Iterative  procedures  but  also  may  produce  an  awkward  result  in 
which  the  estimated  location  parameter  e  is  larger  than  the  smallest  obser¬ 
vation  X^j. 

In  the  present  study,  therefore,  a  convenient  curve  fitting  procedure 
described  below  Is  used  in  conjunction  with  the  least  square  method  under 
the  assumption  that  the  location  parameter  is  known. 


Transforming  Eq.  64  into  the  form 

£n£n{l/(l  -  F^x)}  =  n(£n(x  -  e)  -  £nX) 
and  setting 


y  =  $Mn{l/(l  -  Fx(x)}  v 
u  =  £n(x  -  e) 

X*  =  -nJlnX 

we  can  reduce  Eq.  64  into 
y  =  nu  +  X* 


(65) 


(66) 


(67) 
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which  is  a  linear  relationship  between  y  and  u.  Using  the  observations  x^ 
and  corresponding  cumulative  probabilities  F(x^)  =  i / ( 1  +  n)  listed  in 
Table  1,  we  can  compute  u.  and  y^  which  are  the  values  of  u  and  y  in  Eq.  66 
with  x  replaced  by  x^.j.  These  u.  and  y..  are  used  to  construct  the  square 
error  E2  of  the  form 

E2  =  I  (y-  -  (nu.  +  X*)}2  (68) 

1=1  1  1 

/s  /v 

The  parameter  values  q  and  X*  which  produce  the  least  value  of  E2  are  then 
chosen  as  our  estimates.  Explicit  expressions  for  these  estimates  can  be 
obtained  by  solving  the  equations  3E2/3q  =  3E2/9X*  =  0  as 

n  =  (n  l  u.y,  -  (  l  u,)(  l  y.)}/{n  l  u?  -  (  l  u.)2}  (69) 
1=1  1  1  i=l  1  1=1  1  1=1  1  1=1  1 


X* 


('(  l  u2)(  l  y.)  -  (  l  u.)(  l  u .y. )}/{n  l  u| 
1=1  1  1=1  1  i=l  1  i=l  1  1  1=1  1 


The  scale  parameter  X  is  then  estimated  from  Eq.  66  as 
X  =  exp  (-A*/n) 


(  l  u.)2} 

1=1  1 

(70) 


(71) 


Theresultsof  such  estimations  are  presented  in  Tables  13,  14  and  15 
respectively  for  the  samples  XQPF,  XWPF  and  WPF.  In  deriving  these  results, 
we  have  assumed  that  the  location  parameter  e  =  0,  0.2x^,  0.4x^,  0.6x^ 
and  0.8x^j.  Figs.  19,  20  and  21  plot  y^  against  u^  (although  the  probabil¬ 
ity  and  the  logarithmic  scales  are  used)  respectively  for  XQPF,  XWPF  and 
WPF.  Each  data  point  in  these  figures  represents  x^  -  e  along  the  abscissa 
and  1/(1  +  n)  along  the  ordinate.  Therefore,  at  each  probability  level 
i / ( 1  +  n),  we  see  five  points  corresponding  to  five  different  values  of  e. 
For  each  assumed  value  of  e,  we  evaluate  the  value  of  nio2  from  Eq.  6  using 
the  corresponding  estimates  n  and  X.  In  Tables  13,  14, and  15,  we  also 
list  the  values  of  least  square  error  as  well  as  the  values  of  nw2.  Both 
of  these  values  decrease  as  e  increases.  This  suggests  that  the  goodness- 
of-fit  is  more  satisfactory  if  we  assume  larger  e  values.  However,  it  is 
apparent  from  Fig.  2  as  well  as  Table  2  that  all  the  Weibull  distributions 
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Table  13  Statistical  Analysis  of  Data 
_  JfQPF  for  Welbull  Distribution 


ANALYSIS  OF  EQUIVALENT  INITIAL  FLAW  SIZE  DISTRIBUTION 
IN  THE  CASE  OF  WE  I  BULL  DISTRIBUTION 


FOR  DATA  SERIES  XGPF 


WEIBULL  SHAPE  s  .76818F+00 
wE I  BULL  SCALE  =  ,AA256EfOO 
WEIBULL  LOCATION  =  .OOOOOEfOO 
LEAST  SQUARE  ERROR  =  .72664E+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2  =  .24760E+00 


WEIBULL  SHAPE  a  .75091E+00 
WEIBULL  SCALE  =  .«28?«E+00 

WEIBULL  LOCATION  =  .52000E-02 

LEAST  SQUARE  ERROR  s  .6A314E+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2=  .23008E+00 


WEIBULL  SHAPE  =  ,73159E*00 

WEIBULL  SCALE  =  ,«13l<iE  +  00 

WEIBULL  LOCATION  =  .lOAOOE-Ol 
LEAST  SQUARE  ERROR  a  .58844Et01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2  =  .21037E+00 


WEIBULL  SHAPE  =  .70B78E+00 

WEIBULL  SCALE  s  .39705E+00 
WEIBULL  LOCATION  a  .15600E-01 
LEAST  SQUARE  ERROR  =  .49732E+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2s  .18742E+00 

i 


WEIBULL  SHAPE  s  .67785E+00 
WEIBULL  SCALE  =  .37988E+00 

WEIBULL  LOCATION  =  .20800E-01 

LEAST  SQUARE  ERROR  =  .37987F+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWNPs  .15869E+00 
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Table  14  Statistical  Analysis  of  Data 
XWPF  for  Wei  bull  Distribution 


ANALYSIS  OF  EQUIVALENT  INITIAL  FLAW  SIZE  DISTRIBUTION 
IN  THE  CASE  OF  WE  I  BULL  DISTRIBUTION 


FOR  DATA  SERIES  XWPF 


WE  I  BULL  SHAPE  =  .17920E+01 

WE  I  BULL  SCALE  =  .43272E400 

WE  I  BULL  LOCATION  =  .00000E  +  00 

LEAST  SQUARE  ERROR  =  .50286E+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2  =  .23084E+00 


WEIBULL  SHAPE  =  .16R45E+0! 

WE  I  BULL  SCALE  s  .40981E+00 
WEIBULL  LOCATION  =  .18600E-01 

LEAST  SQUARE  ERROR  =  .4«298E+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2  =  .20895E+00 


WEIBULL  SHAPE  =  .15P9JE+01 

WEIBULL  SCALE  =  .386>53Et00 

WEIBULL  LOCATION  =  .37200E-01 

LEAST  SQUARE  ERROR  =  .37610E+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NwN2=  .1P437E+00 


WEIBULL  SHAPE  s  .14695E+01 
WEIBULL  SCALE  =  .38297F+00 

WEIBULL  LOCATION  s  .5580PE-01 
LEAST  SQUARE  ERROR  c  .30437E+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2s  ,15P92E+00 


WEIBULL  SHAPE  =  .13150E+01 

WEIBULL  SCALE  s  .39009E+00 
WEIBULL  LOCATION  =  .74400E-01 

LEAST  SQUARE  ERROR  s  ,2<l880E  +  0l 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2=  .12513E+00 
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Table  15  Statistical  Analysis  of  Data 
WPF  for  Welbull  Distribution 


ANALYSIS  OF  EQUIVALENT  INITIAL  FLAW  SIZE  DISTRIBUTION 
IN  THE  CASE  OF  WEI0ULL  DISTRIBUTION 


FOR  DATA  SERIES  WPF 


WEIBULL  SHAPE  a  .16945E+01 
WEIRULL  SCALE  s  .87823E+00 
WEIBULL  LOCATION  *  .OOOOOE+OO 
LEAST  SQUARE  ERROR  a  .57231E+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN?s  .30361E+00 


WEIBULL  SHAPE  a  • 1 620RE+0 1 
WEIBULL  SCALE  a  .8fl3«5E+00 
WEIBULL  LOCATION  *  .28000E-01 

LEAST  SQUARt  ERROR  r  .51375E+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2s  .2B299E+00 


WEIBULL  SHAPE  a  .15373E+01 
WEIBULL  SCALE  a  .8088UE*00 
WEIBULL  LOCATION  a  .56000E-01 
LEAST  SQUARE  ERROR  a  .45602E+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2=  .26271E+00 


WEIBULL  SHAPE  a  .14362Et01 

WEIBULL  SCALE  a  .77523E+00 

WEIBULL  LOCATION  a  ,8«000E-01 

LEAST  SQUARE  ERROR  a  .OOSlOEfOl 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2a  .24514E+00 


WEIBULL  SHAPE  a  .12923E+01 
WEIBULL  SCALE  a  .7flfe«8E*00 
WEIBULL  LOCATION  a  .11200E+00 
LEAST  SQUARE  ERROR  a  ,«0969E+01 

GOODNESS-OF-FIT  TEST  STATISTICS  NWN2a  .24091E+00 
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Equivalent  Initial  Flaw  Size  x  (mils) 


CD  CO  lAQQOQQOOQ  O  IA  0) 

o«D  (Dcnoor^coin^n  ni  ^ 


{%)  LS-qeqo^d  SA^BLnuino 
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Fig.  20  Weibull  Probability  Plots  for  Data  XWPF 


[LqeqoJd  9ALqe  [niLing 


Fig.  21  Wei  bull  Probability  Plots  for  Data  WPF 


with  estimated  shape  and  scale  parameter  values  and  with  assumed  location 
parameter  values  (including  those  with  e  =  0)  can  be  accepted  as  the  pop¬ 
ulation  distribution  if  the  significance  level  a  (for  rejection)  is  0.10. 
Indeed,  all  but  one  (WPF  with  e  =  0)  can  be  accepted  even  under  a  =  0.15. 
These  results  are  summarized  in  Table  16. 


Table  16  Best  Estimates  (The  Weibull  Distribution) 


Data 

A 

n 

/\ 

X 

nwn* 

E++ 

£ 

XQPF 

.768 

.443 

.248 

7.27 

0** 

.678 

.380 

.159 

3.80 

.0208+ 

XWPF 

1.79 

.433 

.231 

5.03 

0** 

1.32 

.340 

.125 

2.49 

.0744+ 

WPF 

1.69 

.878 

.304 

5.72 

0** 

1.29 

.746 

.241 

4.10 

- - 

.  1 1 2+ 

*  5%  significance  level  =  .461 

10%  significance  level  =  .347 

**  Assumed  minimum  value  of  e 
+  Assumed  maximum  value  of  e 
Least  square  error 


SECTION  V 


FITTING  TO  THE  PEARSON 
DISTRIBUTION  FAMILY 

5.1  The  Pearson  Distribution  Family 

The  probability  distribution  functions  with  the  following  twelve  types 
of  density  belong  to  the  Pearson  family  [5,  6,  7]. 

Type 

I 

II 
III 

IV 

V 

VI 
VII 

VIII 

IX 

X 

XI 


a(/3+B1  + 

x  (83) 

aW 3+17  -  /ap 

Each  of  these  distribution  functions  is  identified  either  as  a  point. 


Density 

Origin 

Range 

f(x)  =  c  xp(l  -  x/a)q 

x  =  0 

0  <  x  < 

a 

(72) 

f(x)  =  c{l  -  (x/a)2}p 

Mean(=Mode) 

v  11 

X 

Vll 

res 

1 

a 

(73) 

f(x)  =  c  xp  exp(-x/a) 

x  =  0 

V 

X 

V  II 

0 

00 

(74) 

f (x )  =  c{l  +  (x/a)2}"p 
x  exp{-b  tan-1(x/a)} 

Mean 

1 

8 

A 

X 

A 

00 

(75) 

+  %ab/(p-l) 

f (x )  =  c  x”p  exp(-a/x) 

x  =  0 

V 

X 

Vll 

0 

00 

(76) 

f(x)  =  c  xq(l  +  x/a)"p 

x  =  0 

V 

X 

vll 

0 

00 

(77) 

f(x)  =  c{l  +  (x/a)2rp 

Mean(=Mode) 

-00  <  x  < 

00 

(78) 

f (x )  =  c (1  +  x/a)"p 

x  =  0 

-a  <  x  < 

0 

(79) 

f(x)  =  c(l  +  x/a)p 

x  =  0 

-a  <x  < 

0 

(80) 

f(x)  =  c  exp(-x/a) 

x  =  0 

0  <  x  < 

00 

(81) 

f(x)  =  c  x"p 

x  =  b 

V 

X 

vll 

jQ 

00 

(82) 
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a  curve  (possibly  a  straight  line)  or  a  region  in  the  6^  and  62  plane  as 
shown  in  Fig.  22  which  is  constructed  by  extending  the  work  presented  in 
[5,  6,  7],  Fig.  22  plots  the  (b^,b2)  points  for  the  data  XQPF,  XWPF  and 
WPF  and  suggests  that  the  Pearson  distribution  function  of  Type  I  is  sup¬ 
posed  to  represent  all  these  data  well. 

Following  the  procedures  suggested  in  [5,  6,  7],  four  parameters  c. 


a,  p  and  q  of  Type  I  written  in  the  form 

f(x)  =  c  xp(a  -  x)q  (84) 

are  estimated  as 

a  =  h{v2r2^  (85) 

=  %{r1  -  2  +  ^(rj  +  2)(81/r2)3s}  (86) 

q 

c  =  a'(p+q+1)r(p  +  q  +  2)/{F(p  +  l)T(q  +  1)}  (87) 

where 

rx  =  6(82  -  -  1 ) / ( 6  +  36j  -  282)  (88) 

r2  =  6^  +  2)z  +  16^  +  1)  (89) 


and  y2,  and  82  are  given  by  Eqs.  1  and  2  while  r(-)  indicates  the  gamma 
function . 

5.2  Results  of  Estimation 

The  parameters  a,  p,  q  and  c  of  the  Type  I  distribution  function  are 
estimated  and  used  in  Eq.  84  resulting  in  the  following  density  functions. 


For  XQPF, 

f (x)  =  0.0358x"°*956(9.65  -  x)°-0503 

(90) 

For  XWPF, 

f ( x )  =  0.351x'°'562(1.44  -  x)1,43 

(91) 

For  WPF  , 

f (x )  =  0.0187x”°’835(5.07  -  x)1,29 

(92) 
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Unfortunately,  the  resulting  values  of  nw*  are  large:  6.64  for  XQPF,  3.02 
for  XWPF  and  6.03  for  WPF.  These  large  values  of  n^  clearly  indicate  that 
the  Pearson  distribution  of  Type  I  cannot  apply  to  either  of  those  data. 
These  density  functions  are  plotted  in  Figs.  23-25  together  with  the  cor¬ 
responding  normalized  histograms  for  comparison. 
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XQPF 


^ISU9Q  ^iL.LqBqoad 


Fig.  23  Pearson  Type  I  Fit  for  XQPF  Data 
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Fig.  25  Pearson  Type  I  Fit  for  WPF  Data 


SECTION  VI 


FITTING  TO  THE  ASYMPTOTIC 
DISTRIBUTIONS  OF  LARGEST  VALUES 

6.1  Asymptotic  Distributions  of  Largest  Values 

The  following  three  asymptotic  distributions  of  largest  values  are 
tested  for  their  goodness-of-fit  with  respect  to  the  observed  EIFS  data. 

(a)  The  First  Asymptote 

-a^x-g^ 

Fx(x)  =  e“e  -o°  <  x  <  °°  (93) 

where  a^(>  0)  and  g^  are  the  parameters  to  be  estimated. 

(b)  The  Second  Asymptote 

Fx(x)  =  exp[-(x/x2)  ]  0  <  x  (94) 

in  which  a2  (>  0)  and  x2  (>  0)  are  the  parameters. 

(c)  The  Third  Asymptote 

r  a3> 

Fx(x)  =  exp  -(—-^)  (95) 

<3  „ 

with  the  parameters  ct3  (>  0),  x3  (>  0)  and  w(>  0). 

For  each  of  these  asymptotic  distribution  functions,  we  can  find  the 
appropriate  transformations  of  the  probability  scale  and  of  the  random  var¬ 
iate  so  that  the  relationship  between  the  transformed  is  linear,  although 
the  upper  bound  w  in  the  case  of  the  third  asymptotic  distribution  must  be 
known. 

6.2  The  First  Asymptote 

Using  a'procedure  similar  to  that  fer:the  Weibull  distribution  fit. 
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transform  Eq.  93  into  the  form 


y  =  -oijX  -  3j  (96) 

where 

y  =  £Mn{l/F^(x) }  (97) 

Using  the  observations  x^  and  the  corresponding  cumulative  probabilities 
F(x ( i ) )  =  i / ( 1  +  n)  as  listed  in  Table  1,  we  can  construct  y^  from  Eq.  97 
as 


y^.j  =  $Mn{(l  +  n)/i }  (98) 

and  construct  the  square  error  E2  as 

£2  =  iI1{y(i)  '  (-“lx(i)  *  el))!  (99) 

/\  /\ 

The  parameter  values  and  6^  which  produce  the  least  value  of  E2  are  cho¬ 
sen  as  our  best  estimates.  Explicit  expressions  for  these  estimates  can  be 
obtained  by  solving  the  equations  3E2/3a1  =  3E2/361  =  0  and  are  given  by 

Oi  =  -fn  l  x.y.  -  (  l  x.)(  l  y.)}/(n  £  x?  -  (  l  x.)2} 

1  i=l  1  1  i=l  1  1=1  1  1=1  1  1=1  1 

(100) 

=  -((  I  x.2)(  l  y.)  -  (  l  x  )(  l  x.y.)}/{n  l  x?  -  (  l  x.)2} 
1  1=1  1  i=l  1  i=l  1  i=l  1  1  1=1  1  1=1  1 

(101) 

in  which  x^  and  y.  are  written  in  place  of  x^  land  y^  for  simplicity  of 
notation . 

The  results  of  the  parameter  estimation  and  test  of  goodness-of-fit 
are  summarized  in  Table  17  which  lists  the  values  of  the  least  square  E2 
and  nw2  as  well  as  the  best  estimates  of  the  two  parameters  for  the  data 
sets  XQPF,  XWPF  and  WPF.  The  results  clearly  indicate  that  neither  of 
the  data  sets  fits  to  the  first  asymptotic  distribution  at  the  signifi¬ 
cance  level  of  5%  or  10%.  Fig.  26  plots  these  data  sets  on  the  Gumbel 
probability  paper. 
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Table  17  Best  Estimates  (The  Asymptotic  Distributions) 
(a)  The  First  Asymptote 


Best  Estimates 


Data 

A 

“l 

A 

% 

HSfl 

'  E2+ 

XQPF 

.642 

.142 

1.48 

21.3 

XWPF 

.415 

-.106 

2.63 

3.41 

WPF 

.147 

-.629 

.697 

10.9 

(b)  The  Second  Asymptote 


Best  Estimates 


Data 

/N 

a2 

- 

*2 

■n 

Ez+ 

XQPF 

.822 

.113 

.0672 

1.44 

XWPF 

1.87 

.240 

.0306 

1.03 

WPF 

1.76 

.469 

.140 

2.13 

(c)  The  Third  Asymptote 


Best  Estimates 


Data 

A 

a3 

A 

X3 

0) 

(MIL) 

n“n* 

E2+ 

XQPF 

16.3 

30.3 

30 

1.62 

23.0 

XWPF 

105.4 

29.8 

30 

.370 

4.33 

WPF 

40.0 

29.6 

30 

.779 

11.7 

5%  significance  level  =  .461 
10%  significance  level  =  .347 

Least  square  error 
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Fig.  26  Fitting  to  First  A 


6.3  The  Second  Asymptote 

For  the  second  asymptotic  distribution,  consider  the  following  trans¬ 
formation  of  the  random  variable  X  into  Y 

Y  =  Jin  X  (102) 

Then,  the  distribution  function  of  Y  can  be  written  as 

FY(y)  =  e'e  (103) 

where 


B2  =  -a2  £n  x2  (104) 

Analytically,  Eq.  103  is  identical  with  Eq.  93  and  therefore  the  same  method 
of  estimation  as  used  for  the  parameters  of  the  first  asymptotic  distribu¬ 
tion  can  be  employed  in  this  case  as  well:  Eqs.  100  and  101  can  be  solved 

A  A  /V 

for  the  best  estimates  a2  and  B2  for  a 2  and  B2  if  in  these  equations  and 
are  replaced  by  a2  and  B2  and  also  x.  =  x^  by  u^  =  Jin  x^  =  £n  x^. y 

A  A  A 

Once  a2  and  B2  are  found,  the  best  estimate  x2  for  x2  can  be  evaluated 
as 

A  A 

/\  “Bp/Olp 

x2  =  e  c  (105) 

The  estimated  parameters  and  the  values  of  the  least  square  E2  and  nw2 
are  listed  In  Table  17  which  Indicates  that  all  the  data  sets  fit  extremely 
well  to  the  second  asymptotic  distribution:  Observe  the  very  low  values  of 
nw2.  Fig.  27  plots  these  data  sets  on  the  Gumbel  probability  paper. 

6.4  The  Third  Asymptote 
Using  Ja ^transformation  similar  to  Eq.  102; 

Y  =  £n(w  -  X)  (106) 

the  distribution  function  of  Y  can  be  obtained  from  Eq.  95  as 
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86 


ass  ss  §  i!  snn  sn 

d_d  =  (/C)^j  ^i[j.qeqoJd  0AL;B[niun3 

(2a  +  /C2x>)- 
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Fig.  27  Fitting  to  Second  Asymptotic  Distribution  of  Largest  Values 


a3y+e3 

iFy(y)  =  1  -  e'e 

where 


$3  =  -a3  £n  x3 

Transforming  Eq.  107  further  Into  the  following  form 
z  =  a3y  +  33 

with 


z  =  £n  in[l/{l  -  Fy(y)>] 
we  can  construct  the  square  error  E2  as 

E2  =  1I1(Z(1)  -  <V(1)  +b3)}2 

in  which 

z(i)  =  in  £n [(l+n)/(l+n-i )] 
and 


(107) 

(108) 

(109) 

(110) 

(111) 

(112) 


y(^  =  £n(w  -  x^) 


(113) 


A  A 

The  best  estimates  a3  and  83  respectively  for  the  parameters  a3  and  83  can 
then  be  obtained  as  those  values  of  a3  and  $3  that  minimize  the  square  er¬ 
ror  E2  in  Eq.  111.  Hence,  these  best  estimates  can  again  be  obtained  from 


A  A 

Eqs .  100  and  101  by  replacing  oij,  Bj*  x^  and  y^  respectively  by  -a3>  -83, 
yi  =  y(i)  and  zi  =  Z(i)'  best  estl’mate  x3  x3  can  then  be  evaluated 


as 


-83/a3 


(114) 


This  estimation  procedure  requires,  however,  the  knowledge  of  the  up¬ 
per  bound  w.  In  the  present  study,  we  assume  that  the  upper  bound  is  w  = 
l.lX(n),  *‘^x(n)’  3'0x(n)’  Tbe  nun  va^ues  associated  with  these  as- 
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sumed  upper  bounds  are  computed  after  the  corresponding  parameter  values 
are  estimated  and  listed  in  Table  18.  The  table  indicates  that  the  nco2 
values  are  smaller  as  the  upper  bound  w  gets  larger.  Therefore,  for  the 
upper  bound  w,  we  assume  a  crack  size  of  30  mils,  a  value  to  be  used  in 
Section  VII  as  the  (smallest  visible)  crack  size  at  the  time  of  crack  in¬ 
itiation.  The  values  of  the  estimated  parameters,  nw2  and  E2  under  the 
assumption  of  w  =  30  mils  are  then  listed  in  Table  17.  The  table  indi¬ 
cates  that,  for  the  data  sets  XQPF  and  WPF,  the  n«2  values  are  smaller 
under  the  assumption  of  w  =  30  mils  than  those  listed  in  Table  18;  this 
is  consistent  with  the  trend  observed  in  Table  18  that  the  nw2  values  are 
smaller  for  larger  values  of  w.  However,  the  nw2  value  for  the  data  set 
XWPF  under  the  assumption  of  w  =  30  mils  is  larger  than  those  associated 
with  w  =  2.3x^,  ...,  3-0x(n)  in  Table  18  against  the  trend. 

The  following  conclusions  can  be  drawn  from  these  observations:  The 
third  asymptotic  distribution  fits  to  neither  of  the  WQPF  and  WPF  data 
sets,  whether  at  the  5%  or  the  10%  significance  level.  However,  the  dis¬ 
tribution  fits  to  the  XWPF  data  at  the  5%  significance  level  under  the 
assumption  of  w  =  1.6x^  'v  3.0x^  and  w  =  30  mils  while  at  the  10%  sig¬ 
nificance  level  under  w  =  2-7x(n)  ^  3-0x(n)-  How  wel1  the  third  asympto¬ 
tic  distribution  fits  to  the  XWPF  data  under  the  assumption  of  w  >  3.0x^ 
still  remains  to  be  investigated.  No  further  study  has  been  pursued,  in 
this  respect  however,  in  view  of  the  generally  poor  degree  of  goodness- 
of-fit  exhibited  by  the  distribution  for  the  current  data  sets. 

Figs.  28-30  plot  respectively  the  data  sets  XQPF,  XWPF  and  WPF  on 
the  probability  paper  for  Fy(y)  =  1  -  exp[-exp(a3y  +  83)]  when  w  =  l.lx^, 
2.0x^  and  3.0x^.  Note  that  the  abscissa  of  the  probability  paper 
plots  £n(w  -  x) . 
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Table  18  Values  of  nw*  as  a  Function  of  w  (The  Third  Asymptotic  Distribution) 


w 

XQPF 

XWPF 

WPF 

1'1  *  x(n) 

2.28 

0.813 

1.64 

1.2 

2.14 

0.648 

1.43 

1.3 

2.05 

0.566 

1.31 

1.4 

1.98 

0.516 

1.22 

1.5 

1.93 

0.480 

1.16 

1.6 

1.89 

0.454 

1.11 

1.7 

1.86 

0.434 

1.07 

1.8 

1.83 

0.418 

1.04 

1.9 

1.80 

0.404 

1.01 

2.0 

1.78 

0.393 

0.985 

2.1 

1.77 

0.384 

0.966 

2.2 

1.75 

0.375 

0.948 

2.3 

1.73 

0.368 

0.933 

2.4 

1.72 

0.362 

0.920 

2.5 

1.71 

0.356 

0.908 

2.6 

1.70 

0.352 

0.899 

2.7 

1.69 

0.347 

0.887 

2.8 

1.68 

0.343 

0.879 

2.9 

1.67 

0.339 

0.870 

3.0 

1.67 

0.336 

0.864 

=  7.70  mils  for  XQPF,  1.28  mils  for  XWPF  and  3.83  mils  for  WPF 
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Fig.  28  Fitting  XQPF  Data  to  Third  Asymptotic  Distribution  of  Largest  Values 


a_9  -  I  =  (X)Aj  ^.iitqeqoad  9Ai}einumo 
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Fig.  29  Fitting  XWPF  Data  to  Third  Asymptotic  Distribution  of  Largest  Values 
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SECTION  VII 


FITTING  TO  TTCI  COMPATIBLE 
EIFS  DISTRIBUTIONS 

7.1  Compatibility  Between  EIFS  and  TTCI  Distributions 

We  now  follow  Shinozuka  [8]  to  demonstrate  the  existence  of  compatibil¬ 
ity  between  the  EIFS  distribution  function  F^(x)  and  the  TTCI  distribution 
function  FT(t).  The  key  to  the  establishment  of  such  compatibility  is  the 
interpretation  that  the  crack  size  ag  =  a(T)  at  time  T  of  crack  initiation 
is  specifiable  (say,  ag  =  0.03")  and  that  the  crack  size  ag  is  reached  as 
a  result  of  crack  growth  from  the  initial  crack  size  X  =  a(0).  Assuming 
that  the  crack  growth  is  governed  by 

da(t)/dt  =  Q [a (t )] b  (115) 

and  integrating  from  t  =  0  to  t  =  T,  we  obtain  the  following  relationship 
among  X  =  a(0),  aQ  =  a(T)  and  T. 

X  =  a(0)  =  a0/(l  +  aCQ  cQT)1/c  (116) 

where  c  =  b-1  >  0.  With  aQ  specified,  Eq.  116  provides  the  transformation 
necessary  to  establish  the  compatibility.  Indeed,  we  can  derive  on  the 
basis  of  Eq.  116 

Fx(x)  =  1  -  Ft{(x‘c  -  agC)/(cQ)}  (117) 

A  number  of  distribution  functions  have  been  used  for  F^(t).  Partic¬ 
ularly  notable  are  the  log-normal  and  Weibull  distribution  functions.  It 
is  well  documented  that  both  of  these  can  usually  describe  the  observed 
fatigue  data  reasonably  well,  especially  in  the  central  range  of  scattered 
fatigue  data.  They  represent,  however,  widely  differing  underlying  me¬ 
chanical-statistical  models:  The  log-normal  distribution  implies  a  fail¬ 
ure  rate  that  increases  at  first  and  decreases  after  reaching  a  maximum, 
while  the  failure  rate  associated  with  the  Weibull  distribution  is  a  mon- 
otonically  increasing  function  with  time  provided  that  the  shape  parameter 
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is  larger  than  unity.  Since  monotonically  increasing  failure  rates  are 
more  amenable  to  physical  Interpretation,  the  Wei  bull  distribution  is  often 
preferred  to  the  log-normal  distribution.  This  preference  is  sometimes  re¬ 
inforced  further  by  the  fact  that  the  Wei  bull  distribution  has  an  analyti¬ 
cal  form  that  is  mathematically  easy  to  manipulate. 


7.2  Wei  bull  Compatible  Distribution 

Shinozuka  [8]  applied  Eq.  117  to  the  two-parameter  Weibull  (TTCI)  dis¬ 
tribution  and  computed  the  probability  density  and  distribution  functions 

-8 

under  the  assumption  that  Q  =  2.0  x  10"  ,  ag  =  0.04"  and  c  =  4.0.  Yang 
suggests  that  the  same  equation  be  applied  to  the  three-parameter  Weibull 
distribution 

FT(t)  =  1  -  exp{-[(t  -  e)/81a}  t>e  (118) 

in  which  case,  the  Weibull  compatible  distribution  is  obtained  as 


with 


Fx(x)  =  exp{- 


Xjj  =  (a~C  +  cQe)“1/c 


x  <  XU 


(119) 


(120) 


In  the  present  investigation,  the  following  values  of  b  and  Q  are  assumed, 
based  on  some  experimental  evidence. 


Table  19  Values  of  Crack  Propagation  Parameters 


DATA 

XQPF 

XWPF 

WPF 

b 

1.26 

1.26 

1.22 

Q 

2.33  x  10‘3 

2.33  x  10-3 

9.25  x  10"3 

Note:  c  =  b  -  1 
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As  In  the  case  of  the  three-parameter  Welbull  (EIFS)  distribution,  the 
parameter  e  Is  assumed  to  be  known  and  the  shape  and  scale  parameter,  a  and 
3,  of  the  parent  Welbull  distribution  are  estimated  using  the  following  es¬ 
timation  procedure.  First,  transfer  the  random  variable  X  =  a(0)  into 

Y  =  X"c  =  {a(0)}”c  (121) 


then  the  distribution  function  of  Y  can  be  written  in  the  form  of  a  three- 
parameter  Wei  bull  distribution: 


Fy(y)  -  1  -  exp(-(^)°) 

where 


y  >  y 


(122) 


Y  =  xuC  =  aoC  +  cQe 

and 


(123) 


6  =  cQ3  (124) 

Realizations  of  X,  x^j  <  <  ...  <  are  transformed  through  Eq. 

121  into  realizations  of  Y,  y ^  =  x^  >  >  ...  >  y^j  = 

x^.  Note  that  the  largest,  the  second  largest,  ...  of  the  X  sample  be¬ 
come  respectively  the  smallest,  the  second  smallest,  ...  of  the  Y  sample 
as  shown  below: 


Observations  of  EIFS 
Largest  a(0)^ 

Second  largest 

Smallest  a(0)^ 

Defining  Ay  as 

Ay  =  (y^j  -  a"c)/5 


Corresponding  Y 

Smallest  y^  =  {a(0)n>“c 

Second  smallest  y^  =  (a(0)n_1}“c 


Largest 


=  Ca (0)j } 


(125) 


we  assume  that  the  location  parameters  e  of  the  TTCI  (Eq.  118)  and  y  of  the 
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distribution  of  Y  (Eq.  122)  take  respectively  the  following  values 

Ej  =  (j  -  l)Ay/(cQ)  ( j=l >2, ... ,5)  (126) 

and 

Yj  =  3qC  +  (j  -  1)Ay  (j=l »2 . 5)  (127) 

Taking  advantage  of  the  same  procedure  as  used  in  the  case  of  the  Wei  bull 
distribution  fit,  we  Introduce 

z  =  av  +  6*  (128) 

where 

z  =  An  £n{l/ (1  -  Fy(y) )) 

v  =  An  (y  -  y) 

and 

6*  =  -a  An  S  (131) 


(129) 

(130) 


By  means  of  the  least  square  method,  we  can  then  evaluate  the  best  estimates 

A  /V  A  A 

a  and  6*  of  a  and  6  from  Eqs.  69  and  70  by  replacing  n,  A*,  u.  and  y.  there 

A  A  11 

respectively  with  a,  6*,  v^  and  z^  where 

vi  =  £n(y(i)  “  y)  (132) 

and 


z^  =  An  An{(l  +  n)/(l  +  n  -  i)}  (133) 

A  A 

The  best  estimates  8  and  6  of  8  and  6  can  be  found  respectively  from  Eqs. 
124  and  131: 


8  =  6/(cQ) 

A  A 


(134) 

(135) 


Table  20  indicates  the  values  of  nw*  corresponding  to  y.  ( j=l  ,2, . . . ,5) 

n  J 

for  each  of  the  three  data  sets.  In  all  cases,  they  are  smaller  for  smaller 
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Data 


Table  20  Values  of  nw2  as  a  Function  of  Y. 

n 

(Weibull  Compatible  Distribution) 

XQPF 


A  A 


j 

Yj 

nto2 

n 

a 

'  6 

£ 

1 

2.49 

0.869  x  10"1 

2.04 

0.133  x  105 

0 

2 

2.70 

0.941  x  10_1 

3 

2.91 

1.06  x  10"1 

4 

3.12 

1.29  x  10""1 

5 

3.33 

1.89  x  10_1 

Data 

XWPF 

A 

A 

j 

Yj 

no)2 

a 

8 

e 

1 

2.49 

0.345  x  10"1 

4.89 

0.103  x  105 

0 

2 

3.12 

0.381  x  10_1 

3 

3.75 

0.449  x  10"1 

4 

4.39 

0.588  x  10_1 

5 

5.02 

0.950  x  10_1 

Data 

WPF 

A 

A 

j 

Yi 

nu2 

n 

a 

8 

£ 

1 

2.14 

1.51  x  10-1 

4.33 

0.155  x  105 

0 

2 

2.38 

1.61  x  10"1 

3 

2.62 

1.80  x  10-1 

4 

2.86 

2.18  x  10_1 

5 

3.10 

3.14  x  10_1 
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values  of  y  and  hence  for  smaller  values  of  e  indicating  that  the  best  fit 

A  A 

is  obtained  when  e  =  0.  (Table  20  shows  only  the  values  of  a  and  6  for  e 

=  0). 

This  implies,  by  virtue  of  Eq.  120,  that  the  best  fit  is  observed  when 
aQ  (in  this  case  0.03")  Is  taken  as  the  upper  bound  of  the  distribution  of 
the  Initial  crack  size  a(0). 

The  results  of  estimation  are  given  in  Table  21  where  the  best  estimates 
of  a  and  3  and  the  corresponding  values  of  nco*  and  e  are  listed. 


Table  21  Best  Estimates  (Weibull  Compatible  EIFS  Distribution) 


Data 

e 

A 

a 

A 

6 

n“n 

XQPF 

0 

2.04 

1.33  x  104 

.0869 

XWPF 

0 

4.89 

1.03  x  104 

.0345 

WPF 

0 

4.33 

1.55  x  104 

.0151 

Extremely  small  values  of  nw*  above  indicate  that  the  Weibull  compatible 
EIFS  distribution  fits  superbly  well  to  all  the  data  sets. 


Figures  31-33  plot  z^  against  v..  using  the  Weibull  probability  paper 
respectively  for  the  data  sets  XQPF,,  XWPF  and  WPF.  As  in  the  case  of  Figs. 
19-21,  each  data  point  in  these  figures  represents  y.  -  y  along  the  abscis¬ 
sa  and  i/(l+n)  along  the  ordinate  and  therefore  at  each  probability  level 
i/(l+n)  we  see  five  points  corresponding  to  five  different  values  of  y. 


7.3  Log-Normal  Compatible  Distribution 


If  the  TTCI  distribution  F^.(t)  is  a  log-normal  distribution  of  the  form 


FT(t)  =  $ 


fftn(t  -  e)  -  y 


t  >  e 


(136) 
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(%)  A^TfcqBqoad  aATqe-[numD 


83 


Fig.  31  Fitting  XQPF  Data  to  Wei  bull  Compatible  Distribution 
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G 

& 

G 
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Fig,  33  Fitting  WPF  Data  to  Weibull  Compatible  Distribution 


then,  the  corresponding  EIFS  distribution  Fx (x )  can,  by  virtue  of  Eq.  117, 
be  written  as 

fx(x)  =  -^n^x  C  -  .H j  X<XU  (137) 

where  y  and  are  the  quantities  defined  in  Eqs.  123  and  120, respectively. 
Performing  then  the  same  transformation  as  shown  In  Eq.  121,  we  can  derive 
the  distribution  function  of  Y  as 

Fy(y)  =  $[Un(y  -  y)  -  £}/a]  (138) 

where 


£  =  £n(cQ)  +  y 


(139) 


Using  the  same  b  and  Q  values  given  In  Table  19  and  following  the  same 
procedure  as  used  in  the  case  of  the  Weibull  compatible  distribution,  the 
parameters  £  and  a  are  estimated  under  the  assumption  that  the  parameter  e 

A  A  A 

is  known.  Once,  these  estimates  £  and  a  are  found,  the  estimate  y  of  y 
can  be  found  from  Eq.  139  as 


U  =  £  -  in(cQ) 


(140) 


The  location  parameters  e  and  y  are  assumed  to  take  exactly  the  same 
values  as  given  in  Eqs.  126  and  127.  Then,  as  in  the  case  of  the  three- 
parameter  log-normal  distribution  fit  in  Section  III  (see  Eqs.  31  and  32), 

A  A 

we  evaluate  the  estimates  £  and  o  respectively  of  £  and  a  using  the  fol¬ 
lowing  expressions: 

£  =  (1/n)  l  y*  (141) 

1=1  1 


with 


tl  (yf  -  £)2/n]  2  =  [{  l  y*2)/n  -  (£)2] 
i=l  1  1=1  1 


(142) 


-  y) 


(143) 
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Table  22  Indicates  the  values  of  nut-  corresponding  to  y .  (j=l,2,...5) 

*  J 

for  each  of  the  three  data  sets.  As  In  the  case  of  the  Wei  bull  compatible 
distribution,  the  table  indicates  that  (a)  the  nw*  values  are  smaller  for 
smaller  values  of  y  and  hence  for  smaller  values  of  e,  suggesting  that  the 
best  fit  is  obtained  when  e  =  0  and  (b)  the  best  fit  at  e  =  0  implies  that 
the  upper  bound  of  the  EIFS  distribution  is  aQ  (in  this  case  0.03"). 

The  results  of  estimation  are  summarized  in  Table  23  where  the  best 
estimates  of  y  and  a  and  the  corresponding  values  of  nw*  and  e  are  listed. 
This  table  indicates  that  the  goodness-of-fit  is  highly  acceptable;  the 
nw*  values  are  much  smaller  than  0.347,  the  value  associated  with  the  10% 
significance  level,  although  the  Weibull  compatible  distribution  has  re¬ 
sulted  in  even  smaller  noo*  values. 

n 

Figure  34  plots  the  three  data  sets  on  the  Gaussian  probability  paper 
after  they  are  transformed  into  £n(y^j  -  y)  =  £n(x^  -  a^  ) . 

The  upper  and  lower  bounds  Sy  and  of  the  confidence  interval  for 
the  expected  value  £  of  £n(Y  -  y)  in  Eq.  138  can  be  written  as 

=  5  ±  (ti-0#n.i>sMr  (144) 

A  /V 

where  £  =  y  +  &n(cQ),  t^  n  ^  =  two-sided  100(l-a)  percentile  of  the  Student's 
t  distribution  of  (n-1)  degrees -of -freedom  and  s  =  unbiased  standard  devia¬ 
tion.  The  lower  bound  Y^  and  upper  bound  Yy  of  the  corresponding  interval 
for  Y  are  then  obtained  from 

£l  -  £,n(YL  -  y),  £y  =  £n(Yy  -  y)  (145) 

as 

/ 

YL  =  y  +  exp{£  -  (tl  ot  n  l )s/v^n}  (146) 

Yy  =  y  +  exp {l  +  (tl  ajn  jb/vfr}  (147) 
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Tab!  e 

22 

Values  of  nw2  as  a 

Function  of  y. 

(Log-Normal  Compatible  Distribution) 

Data 

XQPF ' 

j 

YJ 

no) 2 

n 

/\ 

V  0 

e 

1 

2.49 

0.184 

9.23  0.55 

0 

2 

2.70 

0.191 

3 

2.91 

0.201 

4 

3.12 

0.218 

5 

3.33 

0.256 

Data 

XWPF 

j 

YJ 

nwz 

n 

/V  A 

y  a 

e 

1 

2.49 

0.109 

9.12  0.23 

0 

2 

3.12 

0.121 

3 

3.75 

0.139 

4 

4.39 

0.168 

5 

5.02 

0.224 

Data 

WPF 

j 

YJ 

no)2 

n 

/N  A 

y  c 

£ 

1 

2.14 

0.181 

9.52  0.26 

0 

2 

2.38 

0.200 

3 

2.62 

0.229 

4 

2.86 

0.275 

5 

3.10 

0.369 
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Table  23  Best  Estimates  (Log-Normal  Compatible  Distribution) 


Data 

e 

y 

o 

"“n 

XQPF 

0 

9.23 

0.55 

0.184 

XWPF 

0 

9.12 

0.23 

0.109 

WPF 

0 

9.52 

0.26 

0.181 

Table  24  The  Lower  and  Upper  Bounds  Xy  and  Xy  (in  mils) 

Corresponding  to  Those  y^  and  yy  for  Several  Confi¬ 
dence  Levels  1-a  (Log-Normal  Compatible  Distribution) 

For  Data  XQPF  (N  =  37) 


1-ALPHA 

XL 

XU 

0.90 

.177 

.340 

0.95 

.161 

.373 

0.99 

.132 

.444 

For  Data  XWPF  (N  =  37) 


1-ALPHA 

XL 

XU 

0.90 

.286 

.373 

0.95 

.275 

.388 

0.99 

.254 

.417 

For  Data  WPF  (N  =  38) 


1 -ALPHA 

XL 

XU 

0.90 

.573 

.762 

0.95 

.549 

.794 

0.99 

.504 

.858 
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Cumulative  Probability;  F  (y)  =  $[{£n(y  -  y)  -  £}/a] 


Log-Normal  Compatible  Distribution 

Fig.  34  Fitting  the  Data  to  Log-Normal  Compatible 

Distribution 


90 


Then,  finally,  the  upper  and  lower  bounds  and  of  the  corresponding 
Interval  for  X  =  a ( 0)  are  evaluated  as 


X 


U 


Y-l/c 
tL  » 


y-1/c 

TU 


(148) 


Table  24  lists  these  upper  and  lower  bounds  for  all  the  data  sets  at  the 
90,  95  and  99%  confidence  levels. 
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SECTION  VIII 


CONCLUSIONS 


Distribution  functions  of  the  Johnson  family,  (including  log-normal 
distribution  functions)  Pearson  family,  Weibull ,  Asymptotic  (the  First, 
Second  and  Third),  and  TTCI  compatible  EIFS  (Weibull  and  log-normal 
compatible  EIFS)  have  been  examined  by  means  of  the  w  method  for  their 
acceptability  in  describing  the.  statistical  characteristics  of  EIFS. 

The  results  are  summarized  in  Table  25  below  where  ’'unacceptable", 
"marginal",  "acceptable"  and  "highly  acceptable"  signify  the  follow¬ 
ing  unless  other  significances  are  indicated  in  the  notes. 


Unacceptable:  nco*  >  0.461;  unacceptable  even  at  the  5%  signifi¬ 
cance  level 

Marginal:  0.461  ^  nio*  >  0.347;  acceptable  at  the  5%  significance 
level  but  unacceptable  at  the  10%  significance  level 

Acceptable:  na^  £  0.347;  acceptable  at  the  10%  significance  le¬ 
vel  . 

Highly  acceptable:  acceptable  at  the  significance  level  higher 

than  20%. 

Notes  (see  Table  25) : 

1.  This  is  the  result  when  a  number  of  values  are  assumed  for  the 
location  parameter  e  and  corresponding  values  of  two  other  par¬ 
ameters  are  estimated.  When  we  treat  all  three  parameters  as 
unknown  and  perform  parameter  estimations,  resulting  nu^  val¬ 
ues  indicate  that  the  Johnson  SL  distribution  is  "unacceptable" 
for  XQPF  and  XWPF.  The  distribution  is  "unacceptable"  for  WPF, 
however,  in  the  sense  that  the  estimated  location  parameter  e 
is  negative  (thus  physically  unacceptable),  although  the  cor¬ 
responding  nw*  value  is  in  the  statistically  "acceptable"  range. 
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2.  This  Is  the  result  when  a  number  of  values  are  assumed  for  the 
lower  bound  e  of  the  distribution  and  corresponding  valiues  of 
the  three  other  parameters  are  estimated.  This  distribution 

is  "unacceptable"  for  XQPF  in  the  sense  that  the  estimated  val¬ 
ues  of  the  upper  bound  X  are  negative  for  all  assumed  values  of 
e. 

3.  This  is  the  result  when  the  location  parameter  e  is  assumed  to 
be  zero.  When  we  assume  that  e  =  0.8x^,  this  distribution  is 
"highly  acceptable"  for  all  the  data  sets. 

4.  Only  Type  I  has  been  tested.  Other  types  are  eliminated  on  the 
basis  of  poor  comparison  between  (b^,b2)  and  (82,82) 

5.  Use  of  these  distributionsrequires  the  knowledge  of  parameters 
b  and  Q  in  the  crack  growth  model  da/dt  =  Qa*3. 


Table  25  Summary  of  Goodness-of-Fit  Tests 


Data 

XQPF 

XWPF 

WPF 

Notes 

unacceptabl e 

marginal 

acceptable 

1 

Johnson  Sg 

unacceptabl e 

unacceptabl e 

acceptable 

2 

Wei  bull 

acceptabl e 

acceptabl e 

acceptable 

3 

Pearson 

unacceptable 

unacceptable 

unacceptable 

4 

unacceptable 

unacceptabl e 

unacceptable 

Asymptote  II 

highly  ac¬ 
ceptabl  e 

highly  ac¬ 
ceptable 

highly  ac¬ 
ceptable 

Asymptote  III 

unacceptable 

acceptable 

unacceptabl  e 

Weibul 1 
Compatibl  e 

highly  ac¬ 
ceptable 

highly  ac¬ 
ceptabl  e 

highly  ac¬ 
ceptabl  e 

5 

Log-Normal 

Compatible 

highly  ac¬ 
ceptable 

highly  ac¬ 
ceptable 

highly  ac¬ 
ceptable 

5 
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